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Abstract: 

• In Section 1, wc present a number of classical results concerning the Generalized Gamma 
Convolution ( : GGC) variables, their Wiener-Gamma representations, and relation 
with the Dirichlet processes. 

• To a GGC variable, one may associate a unique Thorin measure. Let G a positive r.v. 
and Tt(G) ^resp. rt(l/G)^ the Generalized Gamma Convolution with Thorin measure 
t-times the law of G (resp. the law of 1/G). In Section 2, we compare the laws of rt{G) 
and Tt{l/G). 

• In Section 3, we present some old and some new examples of GGC variables, among 
which the lengths of excursions of Bessel processes straddling an independent exponen- 
tial time. 
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Introduction 

1. This survey is concerned with the study of a rich and interesting class of 
infinitely divisible laws on IR+ called the generalized gamma convolutions 
(GGC), a class introduced by O. L. Thorin in 1977 (see [4!)]) and then studied 
thoroughly by L. Bondesson [5]; both the lectures notes by Bondesson and the 
book by Steutel and Van Harn [4.s] contain many results on this class of laws. We 
shall also discuss their close connections to a class of random variables known 
as Dirichlet means whose study was initiated by Cifarelli and Regazzini [10; 11]. 

We shall often make, throughout this paper, the common abuse of language 
which consists of talking about a random variable instead of its law; thus, we 
shall use slightly incorrect terms such as GGC variables, and so on . . . 
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In order to introduce the family of GGC variables as naturally as possible, 
let us consider the 3 sets of r.v.'s (or of laws): Q, S, J defined as follows: 

a) J is the family of infinitely divisible r.v.'s (taking values in R+) 

b) S is the set of self-decomposable r.v.'s, valued in 

c) Q is the set of positive GGC variables 

The Laplace transforms of these variables satisfy: 

£;(e-^^) = exp|-aA- (1 - e-^^)i/(da;)| (A > 0) 

and 

a') \i X ^ J , then a > and v[dx) is a Levy measure, i.e.: 

(1 A x)v[dx) < oo 

b') li X <E S, then a > and iy{dx) — ^ h{x), with h decreasing 
c') U X E Q, then X E S and moreover: 

/•oo 

h{x) = / e--V(d2/) 



"'0 

for a (T-finite and positive measure fi; ^i{dy) is called the Thorin measure 
of X. 

Thus, we have: 

g <ZS<ZJ 

and, of course, the inclusions are strict. 

Another definition of a GGC variable is sometimes given in the literature 
as the limit (in law) of sums of independent gamma variables (with different 
parameters) . 

A comparative discussion of different (but, in the end, equivalent) properties 
of elements of 5 is made in Jeanblanc, Pitman, and Yor [32]. Similarly, it is the 
aim of this survey to gather and to compare different properties for Q. This, and 
more generally, the entire survey is motivated by the fact that, recently, some 
new elements of Q were discovered; likewise, there is some recurring interest in 
the Gamma process (see, e.g. the Festschrift Volume for Dilip Madan [55]). 



2. This survey paper consists of three parts: 

• In the first part, the results being discussed are classical; they are about the 
relationships between different families of r.v.'s and/or processes, namely: GGC 
r.v.'s, Wiener-Gamma integrals, Dirichlet means, compound Poisson processes, 
mixing of Gamma variables, Poisson point processes, GGC subordinators, and 
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SO on. These results arc detailed here in order to ease up the reading of this 
paper for probabilists coming from diverse horizons. 

• In the second part, we discuss a notion of duality for the GGC r.v.'s; in 
particular, when one knows the density, or the Laplace transform of a GGC r.v. 
r, this notion of duality allows to compute explicitely the density, or the Laplace 
transform of the "dual GGC variable". We use this principle to compute the 
Thorin measure of a Pareto r.v. and of a power of a gamma r.v. 

• The third part consists essentially in presenting the explicit computations 
of densities and of Laplace transforms of some particular GGC r.v.'s. In the 
main, these r.v.'s originate from the study of the length of excursion which 
straddles an independent exponential time for a recurrent Bessel process [4] . It 
is noteworthy that, in this third part, the notion of duality presented in the 
second part allows to obtain very easily explicit formulae for the density and 
the Laplace transform of a large number of GGC variables. 

Finally, in the Appendix, we describe an interpolation principle between the 
gamma subordinator (74, t > 0) and a family of GGC subordinators. 

1. Classical results on GGC r.v.'s 

A writing convention. 

1. Each time we write an equality in law between r.v.'s and that on one or 
the other side of this equality, several r.v.'s occur, we always assume that these 
r.v.'s are independent, without mentioning it systematically. 

2. It will be convenient, in some instances, to speak of a r.v. instead of its 
law and vice versa. We hope that no confusion will ensue. 

1.1. The Gamma process 

It is a subordinator - i.e. a Levy process with increasing paths and cadlag tra- 
jectories - 

All through this paper, the reference process is the standard Gamma process 
{it, t > 0), which is a subordinator without drift, and with Levy measure 
^ {x > 0). Thus, its Levy-Khintchinc representation writes (see [2]): 

E[e-^^'] = exp|-iy (l-e-^")^e-"| {X,t>0) (1) 

= exp{-tlog(l + A)} = ^^-i^ (2) 

where formula (2) is obtained from (1) and from the elementary FruUani formula 
(see [34], p. 6): 

roo 1 poo poo \ 

/ (i_e-^")— / e-"V(dz)= / log(l + -)j^{dz) (3) 
Jo ^ Jo Jo ^ ' 
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where, in general, v denotes a positive measure on R_|_, which is cr-finite. Here, 
v^dz) = 6i{dz) is the Dirac measure at 1, but formula (3) shall be useful in the 
sequel. We note that FruUani's formula (3) may be easily obtained by observing 
that the two sides of this formula take the value as A = and have the same 
derivative with respect to A. 

For each i > fixed, 74 follows a gamma law with parameter t: 

P{-1t eda)^^ a*-i da {a > 0) (4) 

This process (74, t > 0) enjoys a large number of remarkable properties which 
make it a "worthy companion" of Brownian motion. In particular, Emery and 
Yor [18] establish a parallel between Brownian motion and its bridges on one 
hand, and the Gamma process and its bridges on the other hand. See also 
Vershik, Yor and Tsilevich [5(1] and Yor ["i'")] for a survey of many remarkable 
properties of the gamma process. 



1.2. Wiener- Gamma integrals and GGC variables 

1.2. a Many times throughout this work, we shall use the properties of the 
integrals: 

/•oo 

r{h) := / h{s) djs (5) 
Jo 

where h : R+ — > M+ is a Borel function such that: 

/"OO 

/ log(l + /i(ii))(iw < CX3 (6) 

Under this hypothesis (6), r{h) is finite a.s. (see Proposition 1.1 below). Of 
course, since the trajectories of the process are a.s. increasing, the integral fea- 
tured in (1.5) may be defined in a path- wise manner, as a usual Stieltjes integral. 
We call r{h) a Wiener-Gamma integral, in analogy with the Wiener integrals 
Jq f{u)dBu,f G L'^{M.-^,du) which constitute the Gaussian space generated 
by Brownian motion (-Bu, u > 0). 

Thus, the family of r.v.'s {T{h) with h such that log (l -I- h{u))du < 00^ 
constitutes the analogue for the process (74, t > 0) of the first Wiener chaos for 
a Brownian motion (-B^, u > 0). 

Let us assume for a moment that, in (1.5), the function h is constant on a 
finite number of intervals, i.e.: 

n-l 
1=0 

for a subdivision = sq < si < ■ ■ ■ < Sn < 00. Then: 

i=0 
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Thus, r(/i) is a hnear combination of independent gamma r.v.'s and we obtain: 

_ n-1 

^(e-^r(ft)) = []£;(e-^'''(^=.+i-^=.)) (A>0) 

i=0 

n— 1 -| 

(from (1.2)) 




expl - V(s.+i-s,) / (1 -e-^ — 







(from FruUani's formula (1.3)) 

Jo y Jo 

(after making the change of variable hiX = y) 



= exp(^-^"(l-e-^^)^ e-^>,(dx) 

where is the image of Lebesgue's measure on M.^ by the application s — > 7^;^' 
This latter formula justifies the following definition: 

1.2.b Definition 1.0: Following ([S], p. 29), we say that a positive r.v. F is 
a generalized gamma convolution (GGC) - without translation term - if there 
exists a positive Radon measure fj. on ]0, oo[ such that: 

E[e-'^] = exp{-/ (1-e-^-)-/ e-->(dz)| (7) 
Jo ^ Jo ^ 

= exp{-^"log(l + ^)M(dz)} (8) 
with: / I log x\ fi{dx) < 00 and / 

J]0,1] J\l,' 



< 00 (9) 



The measure /i is called Thorin's measure associated with F. Thus, from the 
Lcvy-Khintchine formula, a GGC r.v. is infinitely divisible. In fact, since its 
Levy density lr{x) = ^ e~^^fi{dz) satisfies: x — > xl{x) is decreasing, then 
F is a self decomposable r.v. (see, e.g. [37]). Such a self-decomposable r.v. F, 
assumed to be non degenerate, admits a density fr such that fr{x) > for 
every a; > (see [4")] p. 404 ). The study of GGC variables was initiated by O. 
Thorin in a series of papers (see for instance [4!)]). 

1.2.C GGC variables and Wiener-Gamma representations. 
The following Proposition is classical. The reader may refer to Lijoi and Regazz- 



Proposition 1.1 The class of positive GGC variables coincides with the class 
of Wiener- Gamma integrals. More precisely: 
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1. IfT(h) — h{u)d'yu7 then: 

i5(e-A?(/0) ^ cxp| _ log + (10) 

where /i^ denotes the image of Lebesgue's measure on M_|- under the application: 
s — > TT-T' other terms: 

h{s) 



oo 



e" — ds= / e-'^'iihidz) {x > 0) (11) 
Jo 

note that, in (11), h may vanish on some measurable set. 

2. Let r denote a GGC r.v. with Thorin measure ^. Let F^{x) :— Jjp /i(dy) 

(a; > 0) and denote by its right continuous inverse, in the sense of the 
composition of functions. 
Then: 

r^'=^r(/i), with;i(u) = — \ — (12) 

Proof of Proposition 1.1 Let T{h) := h{u) d'ju. It is easily obtained, 
by approximation of by Riemann sums, using also the fact that the Levy 
measure of the process (7^, t > 0) equals — e~^, that: 



00 poc 



iS(e-»™) = cxp{-y duj (1 ^e-} 



cxpi 



after making the change of variable xh{u) — y. We observe, from (11) and (12), 
the equivalence of the conditions: 

J log [1 + h{u))du < 00 4^ j log (1 + fih{dx) < 00 
/ \log x\ fill {dx) < 00 and / 



Hh{dx) 
< 00 



Remark 1.2 



1. Formula (13) may be obtained in a slightly different manner: the process 
(7s — 7s- ■= Cs, s > 0) of jumps of the subordinator (7^, t > 0) is a Poisson 
point process whose intensity measure n equals the Levy measure of (74 , t > 0) 
(see [2]): 

n(dx) = - e'^'dx (14) 

X 

Thus, from the exponential formula for Poisson point processes ([4:!], p. 476): 



E 



exp-A V /(s,e,)] = expj - f ds f {1 - e~^f'^''''^)n{du)\ 
^ ^ Jo Jo 



0<s<t 
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one gets: 

/>oo 

£;[exp-A / h{s)d'-fs] = E{exp~X S h{s){Ajs)) 

Jo s>0 

= exp|-y ds J (l-e-^'*(")")-— dwj, 
which agrees with the expression in (13). 

2. Let h, k : — > M+ two Borel functions which satisfy (6) and assume 

that T{h) ''^'r(A:). Relation (13) and the uniqueness of the Levy measure in 

the Levy-Khintchine representation imply: the images by h and k of Lebesgue's 
measure on are identical. Thus, choosing for k the increasing rearrangement 
h* , (resp.: the decreasing rearrangement ft,*) of h, we obtain that there exists 
essentially a unique increasing function h* , (resp. a unique decreasing function 
h^) such that: 

r(/i)^'=^r(/i*) = r(ft,) (15) 

We recall that the function h* (resp. the function h^) is the unique (equivalence 
class of) increasing (resp. decreasing) function such that for every a > 0: 

mcas (x ; h{x) < a) = meas (x, h*{x) < a) = meas (a; ; h^{x) < a) 

where meas indicates Lebesgue's measure on R+. 

From Proposition 1.1, r{h) is a GGC r.v. and we shall denote by fih the 
Thorin measure associated with T{h). 



1.3. m- Wiener- Gamma integrals, {m,G) GGC r.v. 

In this work, we shall often consider a GGC r.v. whose associated Thorin mea- 
sure has finite total mass. Thus, we shall now particularize Proposition 1.1 in 
this case. 

1.3. a Let to > and h : [0, m] M.^ a Borel function such that: 

/ log{l + h{u))du<oo (16) 

We call TO- Wiener integral of h the r.v.: 

r™(ft) := / h{u)dju (17) 
Jo 

Since (7^, < u < to) and (7m — 7(„j_„)_, < u < to) have the same law, we 
deduce from (17), after making the change of variable m — u = v: 

'^m{h)'^^'^j /i(w)d7„ / h{m~u)d'yu (18) 
Jo Jo 
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We note, in relation with (15) above, that if h is increasing, resp. decreasing, 
then u —>■ h{m — u) is decreasing, resp. increasing. 

1.3.b Let 771 > and G be a positive r.v. such that: 

£;(log+(l/G)) < oo (19) 
We say that a positive r.v. F is a (to, G) GGC if: 

^(e-^^) = exp{ - mj^il - e-^-)^ S(e-"«)} (20) 

Of course, from (7), a (m, G) GGC r.v. is a GGC r.v. whose Thorin measure ^ 
equals: 

H{dx) = mPcidx) (21) 
where Pg denotes the law of G and we have: 

/^(]0,oo[)=m (22) 

Under (21), it is clear that: 

f |logx| ^{dx) < oo and / ^^'^ ^ < oo 

J]0,1] Jll.ool X 

^S(log+(l/G)) < OO (23) 
We denote by Tm{G) the (law of the) r.v. T defined by (20). Hence: 

^[g-Ar,„(G)] ^ exp{-TO^"(l-e-^^)^i?(e-^«)} (24) 

= exp{-TOii;(log(l + (25) 



1.3.C G or 1/G ? How to choose ? 
We have used the notation r,„(G) due to the relation (24). Of course, the relation 
(25) invites, on the contrary, to adopt the notation r„i(l/G). However, we shall 
not adopt this latter notation as the notation Tm{G) is used by L. Bondesson [5] 
who has contributed in an essential manner to the study of the GGC variables. 

1.3.d Proposition 1.1, when the Thorin measure has a finite total mass to, 
admits the following translation. 

Proposition 1.3 r.v. is a {m,G) GGC if and only if it is am- Wiener integral. 
More precisely: 

1) If r™(G)^'^"^f,„(/j), then 

H^) = ^-1,^ , , {ue[0,m]) (26) 
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where Fq denotes the right continuous inverse, in the sense of composition of 
functions, of Fq, the cumulative distribution function ofG. 

2) If r„(/i)*'=^r,„(G), then 

^ ~ Wm) ^'^^ 

where Um denotes a uniform r.v. on [0,m]. 

1.3.e Some classical results 
We gather here some results which are due to L. Bondesson [5] and which we 
shall use in the sequel. Let to > and G satisfy (19). Then, denoting /r„(G) 
the density of TmiG): 

m— 1 

1 (to) 

where g is a completely monotone function (['i], p. 49). Moreover ([■")], p. 50) g 
admits a limit on the right of and: 

5(0+) =cxp{TOS(logG)} (29) 

We note that, since by hypothesis £;(log+(l/G)) < oo, .9(0+) is finite if and 
only if E(\og^{G)) < 00. In Section 2 of this work (see Theorem 2.1) we give 
an explicit form of g when £'(|logG|) < 00. 

• TO may be determined from the knowledge of /p^ (g) '■ 

TO = sup/a > ; lim ■^^"'^^^^'^^ = o] (30) 

(see [5], p. 51). 

1.4- m-Wiener Gamma integrals, m Dirichlet means, Gamma (m) 
mixtures 

1.4. a The preceding discussion leads us to introduce, for every to > 0, the 
Dirichlet process with parameter to, denoted: (Du"', < m < to) and defined 
as: 

(7^1™', 0<u<to)= 0<u<to) (31) 

^ Im ' 

It is well known, and it is an easy consequence of the properties of the "beta- 
gamma algebra" that this process {dIj^\ < u < to) is independent from the 
r.v. 7m, hence from {jy, v > to). Indeed, if 7a and 76 are two independent 
gamma variables with parameter a, resp. b, the basic "beta-gamma algebra" 
states that: 

7a 

,7a + 7fc' 



7a +7b^ {Pa,b,la+b) 
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where Pa,b and 7a+b arc independent and are respeetively beta (a, b) and gamma 
(a + b) distributed. In particular, ^ '^^^ is independent from 7a + 7f,. Thus, for 
every u < m: 

lu _ lu 

Im lu + (7m - lu) 

is independent from 7„ + (7^ - 7«) = 7m 

This allows to write, for h which satisfies (16): 

h{u) di, ^ = ' 7,„ • / ;i(zi) d„ (i^l'") ) (32) 
Jo 

Thus, from Proposition 1.3, we may write for G which satisfies (19): 

r™(G)^'=S,„-A„(G) (33) 

with 

Dm{G) := / — T^— - du(I?l'"^) (34) 
Jo (^^Z™) 

It follows that for every {m,G) GGC r.v., the r.v. Tm{G) is a Gamma (m) 
mixture, i.e. it may be written as: 

r„(G)^'=S„-Z (35) 

where Z is a positive r.v. 

In general, a relationship of the kind: 

X ■ Z^'^^'> X ■ 

(with X and Z independent, and X and Z' independent) does not allow, "via 

simplification" to conclude that Z *■ = ■* Z' . However, when X is a gamma vari- 
able, this "simplification" is licit. More precisely: 

1.4.b The relation (35) determines the law of Z. Indeed, let Z and Z' two 
positive r.v.'s such that: 

(law) , 
1,„- Z = 7™ Z 

Then, for every s € M: 

E[ii^]E[Z'^^]=E[i:^]E[Z"^^] 

Hence: 

E{Z'-') ^ E{Z"') and Z ^ Z' 
1.4.C Remark 1.4 

(We shall not use the present Remark in the sequel of this paper). We come 
back to the notation of point 1 of Remark 1.2 and we denote: (j}™' > J2™'' > 



L.F. James, B. Roynette and M. Yor/GGC and Dirichlet Means 



357 



• • • ) the sequel of the lengths of jumps of the process (7^, u < m) ranked in 
decreasing order. It is not difficult to see that since the intensity measure of 
the Poisson point process ((s, e^), s > O) is ds^-^dx, then the jump times 
{U['^\U^"'\---) constitute a sequence of i.i.d r.v.'s with uniform law on [0,to] 
which is independent from the sequence k> I). Thus: 



h{u)djii 



{law) 



(u/m) 



E- 

A:>1 Fr 



d7„ (from (26)) 

j{m) _ \ ' j{ra) 1 



(36) 



fe>i 



where fc > l) is the sequence of i.i.d r.v.'s with common law 1/G and is 

independent (as a sequence) from the r.v.'s k > I). Indeed: 



P 



-uu" 



< X 



= P 



P 



(rn) 



1 

> - 

X 



>Fg 



= 1 - 



since ^ is uniform on [0, 1]. We deduce from (36) that: 



h{u)du{D^^^) 



1 f™ 

/ h{u) dju ^ 

7m Jo 1^^^ 



(37) 



We note that: '}2k>i ~^ — ~ ^ 



We then define the random Dirichlet measure Pq^^"^ (dx) by the formula: 



J, 



(m) 



k>l 



and we obtain, from (34) and (37): 



xP 



(l/G) 



{dx) 



(38) 



This relation (38) justifies the denomination, for Dm{G), of a Dirichlet means. 

The study of the Dirichlet means can be traced back to an early work of 
Cifarelli and Regazzini [10] which culminates into the more recognized [11]. 
Additional early works on this topic include [14; 21; 26; 56]. 
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See also Bertoin [:]] for an example of (34) where he shows that a Cauchy 
random variable, Ci, may be represented as 

C,^"^'^-- /'cot(7r,s)d7/-"'^ f^M'sds. 



71 Jo 

Where, (A/^, s E [0, 1[) is the right derivative of the convex minorant of a Cauchy 
process. 

1.4.d Multiplication by a beta variable. 
In this section we discuss what happens when D„i{G) is multiplied by certain 
independent beta random variables. This idea, and restatements of the results 
(1.4.d i-iii) below, first appear in James (see [2(S], Theorem 3.1, revised in [-9]). 

1.4.d i) We denote by the set of the laws of r.v.'s of the form 

h{u)du{D^u"^) := D^{h), with / \og{\ + h{u))du<oo, h>Q. 



If Dra{h) € and if Pm^m'-m, with m' > m, is a beta r.v. with parameters 

(m, m' — to) then: 

f3m.,n'-m-Dra{h)eV^'^"> (39) 

Indeed, is independent from 7^' and follows a beta law, with parameters 
(to, to' — to). Hence: 

f3m,m'-m-D„,{h) ^ = ^ h(u) d^^iD^^) 

7m' Jo 

{law) Jrn , / n dju 

— ' n[u) ■ 



7m' Jo 7m 

c; ~i 1 Z"™' 
- — / h{u)l[o,rn]{u)d^iDi"''^)eV("^") 
7m' Jo 

1.4.d ii) In the same spirit as for the preceding point, we note that, if G is 
a positive r.v. such that i?(log^G) < 00 and if to' > to, then: 

where, on the RHS of (40) , G and Yp are independent and Yp is a Bernoulli r.v. 
with parameter P — 

P{Yp^l)^p^l-P{Yp^O) (41) 
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Indeed, wc deduce from (24): 

^(g-Ar„(i/G)) ^ exp|-m^"(l-e-^-)^i?(e-^/«)| 



exp(-m' /"(I - e-^-) — — Eier^/^) 
[Jo X m' 



1.4.d iii) We now write the relation (40) in a slightly different manner with 
the introduction of the r.v.'s D,n{l/G) and £'?ri(g^)- We have, from (40) and 
(35) for m' > to: 

r„(l/G) = 7™I?™(l/G') = ™'Vg1^/ ■ 



Hence: 



lm'D^(l/G) = 7m' Ai 



7m' ^ GYp 

so that, from point 1.4.b: 



Tr, 
7r)^ 



.i?„.(l/G)^'=)A„'(^) i.e. : 

/3m,m'-m ' A„(l/G) ^'^"^ Z?„. (^) = ^) (42) 

In particular, for to < 1 and m' = 1: 

(3^^^_^.DUl/Gf=''DA-^) (43) 

1.4.d iv) Some elements of P^™). 
Let T denote a positive r.v. which belongs to the Bondesson class B (see [ ,], p. 
73, Th. 5.2.2) i.e. whose density /t writes: 

/t(.t) = C xf^-^ hi{x)h2{l/x) (/3 e R, a; > 0) 

with ^ 

hj{x) =e^p\ -bjx+ log( — ■ — )vj{dy)\ .7 = 1,2 
^ Jo \x + y/ J 



and 



< oo 



1 + y 

Then, from Bondesson ([•")], Th. 5.2.2, p. 79), we know that: 

:= 7m • T (44) 
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is GGC (for every m > 0) with some, possibly unknown, associated Thorin 
measure we denote as Assuming furthermore that £^(T^"') < oo, we get, 
from (30): 

/r,„(^) 

r:\0+ X° 

But. an elementary computation, starting from (44), shows that: 



Ai„(]0,cx3[) =sup{a>0; lim %M = o| (45) 



Thus, we have: 

Mm(]0, ooQ = m 

Hence, there exists, from Proposition 1.3, a positive r.v. Gm such that 
£^(log^ (^T— )) < oo and also such that, from (33): 

Tm ^'^Sm • T = r™(G™) = 7„ • D^{Gjn) (47) 
Thus, from point 1.4.b 

T^'=^An(G,„)Gl?(") (48) 

We now summarize what we have just obtained; 
Proposition 1.5 

Let T denote a positive r.v. which belongs to B, such that: 

E^T-m'^ < oo 

Then: 

1) T e P^™^ 

2) For every m' > to, • T <E P^'"') (from (^5)) 

In particular. Proposition 1.5 may be applied in the following cases: 

• If r is a generalized inverse Gaussian r.v., i.e. its density is given by: 

frix) = Cx'^-^expl - ^(^12^ + ^)} ' l[o,oo[(2;) 

(/3 € M, ci,C2 > 0) then, for every m > 0, T G 2?(") (ahhough T is a 
GGC variable with Thorin measure of infinite total mass (see [5], p. 59). 

• If T is a Gamma r.v. 79 with parameter 9 then for every to > 0, jo & P^™) . 
Indeed 751 S B and for every in < 6: E(-^) < 00. On the other hand, 
assuming still m < 9, we have: 

An,e-™ • le ' Ira G V^'>^ from (40) 
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• If r is a positive stable r.v. with index a: 

S(e-^^) = cxp{-CA"} (A,C > 0,aG]0,l[) 

then, for a < i, T G 1 (sec [r.], p. 85-88) and hence T e X>^™) for every 
m > since < oo. 

Remark 1.5 Epifani, Gughclmi and MehUi ([16], section 4; see also [17]), posed 
the natural question of which kind of probability measures are the laws of Dirich- 
let means. They were able to find some examples in cases where those particular 
random variables possessed all finite moments. One sees that Proposition 1.5, 
in a rather simple way, identifies a large number of possible distributions. 

1.4.e Another representation of TmiG). 
We have been interested mainly until now in the distributions of the (m, G) 
GGC r.v.'s. We shall now describe a realization of such a r.v. with the help of 
a compound Poisson process. Besides, this realization allows us to show that a 
(m, G) GGC solves an "affine equation" . For a nice survey of these equations, 
see Vervaat [51]. 

Let m > and let X be a positive r.v. We shall say that {Yt, t > 0) is a 
(m, K) M-|_ valued compound Poisson process if: 

i=l 

where (i^i, K2, ■ ■ ■ ) is a sequence of i.i.d. variables, distributed as K , and with 
{Nt, t > 0) a Poisson process with parameter m, independent of the sequence 
(A'i, i = 1, 2 • • • ). In particular, Nt is a Poisson r.v. with parameter mt. 

Proposition 1.6. 

Let r™(G) a (to, G) GGC r.v. {with E{\og+{^)) < 00 and m > . Define K 
by: 

X *■ = ■* — (e, a standard exponential r.v. independent ofG) 

Then 

1) r,„(G)^'"="'^ He-'dYt 
Jo 

where {Yt, t > 0) is a {m, K) compound Poisson process. 
2) r„j(G) satisfies the affine equation: 

r„(G)''=^[/i/™[r„,(G) + i^] 

{^with U uniformly distributed on [0, 1]). 
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Proof of Proposition 1.6. 

i) We first prove point 1. We consider (Yt, t > 0) a. {m,K) compound Poisson 
process. Then, approximating e~*'dYt by the Ricmann sums 

we obtain: 

E (^exp{ - xj\-'dYt]^ = cxp{ - ^"(1 - e~^''V([«,oo])^} 

where ^ is the Levy measure of the subordinator (Yt, i > 0). Since this subor- 
dinator is a (to, K) compound Poisson process, we have: 

/i([t;,oo[) = m P{K > v) ^ m P (^-^ > 
= mP{e>vG)=mE{e-'"'^). 

Hence: 

£; (^exp| - A^°°e-*drt}^ = exp| - to^°°(1 - e"^^)^ £;(e-"'^)| 

ii) Wc now prove point 2. We have: 



Jo 



e~*dYt = / e~*dYt + / e'^dYt 

Jo JTi 



(where Ti is the first jump time of {Nt, t > 0)) 

/•oo 

= e-^'K,+e~^' / e-'d{YT,+t) 
Jo 

and wc observe that: 

Jo Jo 

and e~*(i(lTi+t) is independent of Ti. 

J. 5. T/ie subordinators {Tt(G), t > O) 

Until now, we have been interested uniquely in the "individual" GGC variables. 
However, to each GGC r.v. F we may, since P is infinitely divisible and posi- 
tive, associate a unique subordinator (Pj, t > 0) such that Pi r. It is this 
subordinator which we shall now define and describe. 
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1.5. a Let T denote a GGC r.v. From Lcvy-Khintchine formula, there exists 
a (miique) subordinator (Ff, t > 0) such that: 

Fi^'=^F (49) 

hence we have: 

E{e-^^') =cxp^-t J (l-e"^")^y e-^^Xdx)} (50) 

where /z denotes the Thorin measure associated with F. 

1.5.b Let G denote a positive r.v. which satisfies (f9). Then, there exists a 
subordinator (Ft(G), t > 0) which is characterized by: 

dx 

X 

A- 



£;(e-^r,(G)) ^ c^p^-tj (l-e-^-)^iJ(e-G)| (51) 
= exp{-ti?(log(f + (52) 



In particular, for every t > 0, Ft(G) is a (i, G) GGC r.v. Thus, there exists, 
following (33), a family of r.v.'s Dt{G), t > 0, whose laws are characterized by: 

Ft(G)=7fA(G) (53) 

and, for every i > 0, from (34) and Proposition 1.3: 



, , (law) /"* 1 , (law) /"* 1 , 



^G'(l-f) 



(54) 



A(G) ''-^ f -1- d..[D^^) ^'-^ f ' d^iD^) (55) 

We note that the relations (54) and (55) are only true for fixed t, for any t > 0, 
but do not hold as equalities in law between processes. On the other hand, since: 
£;(g-A7t-) _ (i+x)t ; deduce from (53) that: 

^(e-Ar.(G)) 1 \ (5g) 

\{l + XD,{G))') 

1.5.C Some elementary properties of (Ft(G), i > 0) and (^Dt{G), t > 
Let G denote a positive r.v. which satisfies (19). 

Proposition 1.7 

1) The family of laws of the r.v. 's Dt{G), t > solves the equation: for every 
t,s> 0: 

"ft+s ■ Xt+s 7t • + -jsXs (57) 
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2) The family of laws of the r.v. 's [Dt{G), t > O) solves the equation: for every 
s,t> 0; 

Xt+,{G) ^ /3t^,Xt + (1 ~ (3t^,)X, (58) 

3) We assume furthermore that E(^^^ < oo (^which implies E (log^(^) < oo) . 
Then: 

ii) for every t>0, E{Tt{G)) ^ t E{l/G) and E{Dt{G)) = E{l/G) (59) 
3m) -rt{G)^E{l/G) and DtiG) ^ E{1/G) (60) 

t t^O t^OD 



3m) [r,(G)]'^3t/ and [AlG)]^^^! (61) 
where U is uniform on [0, 1] . 

Point 2 of this Proposition 1.7 is due to Hjort and Ongaro (sec [27]), which can 
be seen as a consequence of Ethier and Griffiths [19, Lemma 1]. 

We note the following remarkable feature of points 1 and 2 of Proposition 
1.7: the affine equations (where the unknowns are the laws of the (Xt, t > 0): 

7t+s • Xt+s *'='^ jt-Xt+-/s- Xs {sit > 0) 

and Xt+s *'=^ Pt,sXt + {l- f3t+s)Xs is,t > 0) 

both admit infinitely many solutions: Xu = Du{G), u > for every r.v. G 
which satisfy (1.19). 

Proof of Proposition 1.7 

Point 1 : 

7t+. • A+s (G) ^ 7t DtiG)+ 7. Ds (G) (62) 

follows from: 

Tt+siG) ^'^^^Tt{G)+Ts{G) 
since (TtiG), t > O) is a subordinator and since, from the definition of Da{G): 

Ta{G) A(G). 
We now show (58). From the beta-gamma algebra, we have: 

(7t, 7s) ^'='^ {Pt,s 7t+s, (1 - fit,s)lt+s) (63) 
hence, plugging (63) in (62), we obtain: 

-It+s ■ Dt+s{G) ^'=^Pt,s It+s Dt{G) + (1 - f3t,s)lt+sDs{G) 
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which, using point 1.4.b, implies relation (58). We now prove 3i). We have: 



A=0 



d_ 

dX 
tE(l/G) 



^exp{-t^;(log(l + A))}|^^^ 



and we deduce from (53), that: 

E[Tt{G)] = E[jt ■ A(G)] = tE{Dt{G)) =tE{l/G) 

We now prove 3 ii). The a.s. convergence (hence, the convergence in law) of 
i ^t{G) as t ^ oo towards E{\/G) follows from the law of large numbers. 

We prove that Dt{G) — > E{l/G). Indeed, we may write: 



Tt{G) {law) 



t 



D 



:(G)(^-l)+A(G) 



(64) 



and we shall prove that: 



A(G) 



7t 
t 



which implies that A(G) E{l/G). Since E{Dt{G)) ^ E(^), this implies 

that the family of laws of {Dt{G), i > O} is tight, hence convergence in proba- 
bility implies convergence in L^. It remains to prove: E (Dt{G) 1^ — l|) — > 0. 
But we have: 



E[Dt{G) 



t 



^El-)E 



t 



<El-]lEl {^] -1 



It 



G 



t t^c 



Hence the result. 

We now prove point 3 in). 



If we knew that (Z)t(G))* > 1, then, combining this result with the classical 



one: (74)* *— U, we would deduce: 







{r,{G)y^'=>{^,Y-{D,{G)y'^>u 

In fact, we shall proceed in the other direction. We shall show further (see point 
7 of Remark 2.2) that (rt(G)) ^'^^U. The relation (rt(G))* (7*)* (A(G))* 
and: 

(7t)* U, (rt(G))* U then imply easily that: 
{D,iG)y'^h. 



, t {lau 
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Remark 1.7 Lijoi and Rcggazini [M>] have shown that the support of the law 
of Dt{G) is the closure of the convex hull of the support of the law of 1/G. In 
particular: 

• EG <a a.s., then A(G) > ^ a.s. (65) 

• If G > a a.s., then A(G) < - a.s. (66) 

a 

In Section 3 of this survey, we shall verify this assertion on inspection of numer- 
ous examples. 

1.6. Some examples of GGC subordinators 

Let r denote a GGC variable with associated Thorin measure fi, and let (Ft, t > 
0) denote the subordinator such that Fi = "* F. We have: 

r r°° rl-r 1"°° ^ 

Eie-^^') :=^cxp\^-t J (l-e-^^)^y e-^Xdz)| (67) 

Such a subordinator is called a GGC subordinator. 

Here are now some examples of such subordinators. They are lifted from a 
paper by H. Matsumoto, L. Nguyen and M. Yor [.'s] on one hand, and from the 
study of hyperbolic subordinators made in Pitman and Yor [i ] on the other 
hand. The reader may refer to these papers for further information. In the study 
of these examples, we shall denote subordinators with curly letters, especially 
to avoid some possible confusion with the modified Bessel functions, which are 
traditionally written with ordinary capital letters. 

1.6. a The hyperbolic subordinators (see [41] for a probabilistic description 
of these subordinators.) 

i) The subordinator [Ct, t > 0) is characterized by: 

E{e-^'^' ) = cxp{ - t log cosh V2A} = ( ^-^) * (A, t > 0) (68) 

^cosh v2A^ 

Its Levy density Ic equals: 

, oo 2 

Hx)^-Y,c^p(^~^i2n-lfx) ix>0) (69) 

71—1 

Hence, its associated Thorin measure, (i.e.: the Thorin measure of Ci) 
equals: 

oo 

fic{dx)^J2^4i2n-lAdx) (70) 
n=l 

It has infinite total mass. 
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a) The subordinator [St, t > 0) is characterized by: 
Its Levy density equals: 

-.00 2 

n— 1 

Hence, its Thorin measure equals: 

00 

fisidx) = V (dec) (x > 0) (73) 

n=l 

It has infinite total mass. 

We note that the subordinator {%, t > 0) which is characterized by: 

= ( . f^Jl^y = /tanW2Ay 

VcoshV2A/ V V2A / V V2A / 

satisfies: 

{Cu t>Q)^'"^\St+Tt, t>0) (75) 

and that its Levy density equals: 

It{x) = lc{x) - ls{x) 

However, this subordinator t > 0) is not GGC, as its would be 'Thorin 
measure' fix is a signed measure: 

1.6.b The subordinators (J^'"', t > 0) and (/cf \ t > 0) 
We denote by 1^ and i^i, the modified Bessel functions with index v (see [34], 
p. 108). 

i) The subordinator {Ji^\ t > 0) is characterized by: 

E{e~^^'') = (1 + A+ v/(l + A)2 - 1)-* = exp{-t ^ {l-e'^'^)^ I,{x)e-^] 

(76) 

Its Levy density Z j(o) equals: 

;j(o,(a;) =/o(x) ^ (a;>0) (77) 

and its Thorin measure Hjio), with total mass equal to 1, equals: 

1 dx 

fij(o){dx) = , . , l[o,2] {x) (78) 

TT y^X(2-X) 
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This formula was obtained by Bondesson ([7], Ex. 5.2). In Section 3, Theorem 
3.1, we shall meet again this subordinator {Ji^\ t>Q). 

ii) The subordinator {icf \ t > 0) is characterized by: 

E{e-^'^^°^) = exp| - ^ (argcosh(l + A))^| 

= exp|-ty (l-e-^")^ifo(x)e-"| 



(79) 



Hence, its Levy density equals: 

/yc(o) {x) = Ko{x)- — x>0 (80) 

X 

and its Thorin measure Hjcm equals: 

dx 

fMK:io){dx) = l[2,oo[(a;)rfa; (81) 

V x[x — 2) 

We note that this latter formula follows from: 

/•OO 

e-^'Koix) = e-" / e-"(^°"h*)dt (see [34], p. 119) 



e 



-xu _ f —XV d^ 



^v? - 1 Ji ^v{v ~ 2) 



Hi) The subordinators {K^f^ ^ i > 0) and {jf^^ i > 0) are connected via the 
subordination relation: 

(:7(°',t>0)''?.-'(4";^^„,t>0) (82) 

where, on the RHS from (82), the processes (/ci°\ u > 0) and (51/2 (i), * > 0) 
are independent, and where (S\i2^\ i > 0) is the stable subordinator with 
index 1/2 characterized by: 

£:(e->Si/2(t)) = cxp( - i\/2A) (83) 

iv) The v.a. /c|'^\ for t fixed, may be realized in the following manner: 

let (6„(s) ; < s < u) the brownian bridge with length u (with 6„(0) — bu{u) = 

O) and let: 

A{bu) := / exp(25„(s))ds 
Jo 

Then: 

^""^ [Aih/.)] ^ t[[ cxp (^) ds) (84) 
where, in (84), (6(s), < s < l) denotes the standard brownian bridge bi. 
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v) In relation with the preceding. D. Dufresne and M. Yor [1^], in a work in 
preparation [15], establish the following formula: let (blf\s), < s < u) denote 

the brownian bridge with length u, starting from and such that b\^\u) = x. 
Let: 

A[bi^^ := / exp{2bi^Hs))ds (85) 
Jo 

Then, D. Dufresne and M. Yor obtain the representation as a Wiener-Gamma 
integral of the r.v. : 



1 {law) 



A[bl^''] Jo e^"' + 2e^ cosh {su) + 1 

and in particular: 



djs (86) 



1 {law) f°° djs ^g^^ 



A[bu] Jo 1 + cosh (su) 

1.6.C The subordinators (J'/''\ t > 0) and {ic[''\ t > 0) 
They are obtained by replacing in (76), resp. (79), Iq by I^, resp. Kq by K^. 
i) The subordinator {Jt^\ t > 0), which is defined for v > —1, is characterized 
by: 

E{e-^^i ) = expj ~t j (1 - e-^") -- e-"/,(a;)| (88) 
Thus, its Levy density l-jiv) equals: 

lJ^^■^ {x) = ^I^{x) x>0 (89) 

Its Thorin measure, in the case —1/2 < i/ < 0, is obtained by the following 
calculus: 

^•^^'^ = r-J}^\ 1 m t^r^'/^e^' + e-^')dt (Re,. > -1/2) 

yhTT{iy+ 1/2) Jq 

(see [34], p. 119, form.5. 10.22) 

and: ^ 

e-^/^(z) = az'' / (l-i2)''-i/2(g-^(i-*) +e-^(i+*))dt 



But, for 1/ < 0: 

^ /-OO 



r(-^^) Jo 



s 



-u-l^-z{a+s)^^ 
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Hence: 



I,{z) = Cl s-'^-i(e--'(i-*+^) +e---(i+*+^')ds 



e-'''{h - 1 - t)-"-'dh 



i+t 



C' 



dh / dtil^t^Y-^'^lih-l+ty-Hhyl-t 



+ {h-l- t)- 



Mi>l+t I 



Hence, in the case — 1/2 < < 0, the Thorin measure of {J^''\ t > 0) admits 
a density equal to: 



fij,.,{dx)/dx^ cu {i-t^r-'^mx-{i-t)) " 'i,>i^t 

Jo 

+ {x - {1 + t)y^ha,>i+t}dt 

Hence, if -1/2 < < 0, {Jl", t > 0) is a GGC subordinator. 

a) The subordinator {ICl'^\ t > 0), defined for i/ < 1, is characterized by: 

Its Levy density I^^m equals: 

iKi^^ix) ^ — K,{x) {x>0) 

X 

Its Thorin measure /i^c*"' ' obtained by using the formula {[:]-\], p. 119) 
K^{z)= e"^'=°"^"(coshi^u)dw 



(90) 



(91) 



equals: 



{dx) ~ cosli[i^arg cosh(a; — 1)] 



dx 



Vx{x ~ 2) 



l[2,oo[(a;) (92) 



The formula in (92) simplifies using, 



cosh[i/argcosh(x-l)] = -^{{{x - 1) + \/{x - 2)x) +{{x - 1) + \/{x - 2)x) }, 



smce: 



arg cosh(y) = log (y + - l) . 
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1.7. The subordinator {VtiSp), t > O) (0 < (3 < 1) 

We give here a realization for the subordinator (TtiSfs), t > O), where is a 
positive /3-stable r.v. 

Let a > 0. We denote at^n '■= {si = < S2 • • • < s„ = t} a subdivision 
of the interval [0,i\ with mesh S{at^n) ■= sup (s^+i — s,). The limit along the 

i 

decreasing filtering set of sequences crt.„ whose meshes d{at,n) tend to zero, of 
Es,e^t,„(T«fc+i - ^s*.)" '^'^ists a.s. Observe that: Es,e^,,„ (Ts^+i - 7sJ") = 
E, P. r(a+..+i-.,) ^ ^p^^^ Let: 

:= lim ^ (7.,+, - = E " T-)" (^3) 

o ( ''"t J7 ) — 'U 

Sfc6'7t,n 0<s<t 

Let 

From Section 1.4, A^t^^; is a Poisson r.v. with parameter t ■ /T/o - — du. But, 
evidently: 

iV*,, = #{s<i; (7s-7s-)" >4 (94) 
We deduce from (94) that the Poisson point process 

I (7s -7s-)", s > 0) admits n^°'^{dx) = - — e~^'^° (x > 0) 

ax 

as intensity measure. Hence (^'■"■'(i), t > O) is a subordinator whose Levy 
measure equals n^"\dx): 



ij(e-Ay<°'(t))^cxp 



I a Jo X J 



This subordinator, for every a > 0, is self-decomposable but it is GGC only in 
the case a> 1. In this case: a > 1 we have: 

£;(e-^^'°H*)) = cxpf - - - e-^-)— £;(e--^V=)\ 

\ a Jn X J 



where Si/a is a positive ^ stable r.v. In other terms, with the notation of Section 
1.5.b: 

(rt(^;3), i>o)^'^"'(y(")(to), t>o) (« = ^>i) (95) 



1 . 8. Remark 



We now end up this Introduction by indicating how the study of the GGC 
subordinators may be embedded in a more general one. 
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A subordinator {Nt, t > 0) is said to belong to the Thorin class T^^^S.^), 
with X > (see [Grig]), if its Levy measure admits a density In of the form: 

iNix) = x^-'^k{x) x>0 (96) 
where fc is a completely monotonic function, i.e. it may be represented as: 

k{x) := / e-^y fi{dy) (97) 

for a positive Radon measure ^ carried by M!j_. Thus, the subordinators which 
we study in this work, i.e.: the GGC subordinators, belong to the class T(^)(M+). 
The class r(2)(IR+) has been studied by Goldie [24], Steutel [47] and Bondes- 
son [G]. The r.v's which belong to this class are the generalized convolutions of 
mixtures of exponential laws. Note that T^^\R+) C T^^">{R+) if x < x'- We 
also note that [BNMS] present extensions of these notions to M''. 

In the same manner as condition (9) is necessary and sufficient for a measure ^ 
to be the Thorin measure associated to a positive r.v., B. Grigclionis [2-'i] obtains 
a necessary and sufficient analytical condition so that a measure ^ defines, via 
(96) and (97), a subordinator {Nt, t > 0) which belongs to the class T(x)(R+). 

1.9 We now detail the contents of the sequel of this paper: 

• In Section 2, we present a duality result which connects on one hand the r.v.'s 
Tt{G) and A(G) to the r.v.'s Tt{l/G) and Tt{G) on the other hand. 

• In Section 3, we study in depth the examples of subordinators Tt{Gra) (0 < a < 
1) for which we know how to compute explicitly their Laplace transforms, i.e.: 
their Levy exponents, as well as their densities at any time, and their Wiener- 
Gamma representations. 

2. A duality principle 

Throughout this section, G denotes a positive r.v. such that i?(|logG|) < oo. 
Thus, we have: 

£;(log+(l/G)) < C30 and E {log+ (G)) < oo (98) 

Consequently, the subordinators (Tt{G), t > O) and (ri(l/G), t > O) are 
well defined. We denote by tpa (resp. V'i/g) the characteristic exponent (i.e.: 
the Bernstein function or Levy exponent) of the subordinator (r((G), t > O), 
resp. (rt(l/G), t > 0). 

^g(A) := / (1 - e-^n — i?(e--«) (99) 
Jo a; 

and the same formula for V'i/g obtained when replacing G by 1/G. We denote 
by Fg the cumulative distribution function of G and by F^^ its right continuous 
inverse, in the sense of the composition of functions. 



L.F. James, B. Roynette and M. Yor/GGC and Dirichlet Means 



373 



2.1. The duality theorem 
Theorem 2.1 (Duality) 

1) F-/^{y)-F^\l-y) = l a.e. {y & [0,1]) (100) 

2) For any A > 0: 

^i/g(A)-Vg(1/A) = £;(logG)+logA (101) 

3) /rt(i/G) the density of the r.v. Vt{l/G), equals: 



/r.(i/G)(-)=e-*^('°s«)i?{(^iM)'^'j,_,(2y^r;(G))} (102) 



'rt(G)' 

X 

where J^, denotes the Bessel function with index v: 



oo 



'^-(-)-EM) %(fc + i)r(fc + . + i) y ' N<-'|A>-g^l<- 

(103) 

(see [34], p. 102) 

»*)/r.(i/G)(a:)-^e-*^(iogG)^(g-xD.(G)) > q) (^04) 



^^e[ ^ Ue*^(i°gG) (105) 

V(A(G))V 

^) '^^{t^ ^77^71)= cxp{-i(V'i/G(A)-i?(logG)} (106) 

\(A + l?t(G)) / 

a) The densities Jd^^q-j O'l^'d fotii/G) of, resp., Df{G) and Df{l/G) are 
related by: 

fn,(i/G){x) = x'-^ e-*^(i°s«) fD,(G){llx) [x > 0) (107) 
Remark 2.2. 

1) The relation (101) is found, under a slightly different form, in Bondesson ([5], 
p. 48, 111 "Curious composition"), under the only hypothesis iiJ(log'''(l/G)) < 
00. We have chosen to work under the stronger hypothesis i?(|log G|) < 00 since 
it seems that only under this hypothesis can we obtain the most interesting 
results. 

2) Of course, all formulae of Theorem 2.1 arc "involutive" , in that G may be 
replaced there by 1/G. 

3) Formula (104) agrees with, and makes more precise, the result of Bondesson 
which we recalled in (28) and (30). In particular: 
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4) Formula (101) may be generalized as follows: 

Let a, 6, c, (i e M**, with ad — be ~ ±1 and let, for a;, A > 0: 

ax + b ^ dX + b 

so that <t{G) be a positive r.v. Then: 

-V'g(2^(A)) + 7/',(G) (A) = log (cA + a) + fc 

with fc := i?(log(^^y^^ • ^^^)) — log (c Ao + a) where Aq is a fixed point of a. 
The relation (101) corresponds to a = d = 0, b = c = 1. We shall study, in the 
Appendix, the case where a ~ d ~ sinhu, 6 = c = coshu (u > 0). 

5) Cifarelli and Regazzini [; 1], Cifarelli and Melilli [0] have obtained the density 
of Dt{G) for t>\ and James, Lijoi and Priinster [M)] have obtained it for t < 1. 
For t < 1, they obtained: 

/d,(G)(.^) = j\x - uY-'^{6,{u))du (110) 

with 

Otiu) := i sin (t TT i^i/G(^i)) exp{ - t E (log - i) 1„^^) } (111) 

The proof of this formula (110) hinges upon the knowledge of the density of the 
r.v. which is defined by: (see (41) and (42) with m = t and m' = l). 

A,i-*-A(l/G)<'='z?,(^) (t<l) 
This density equals: 

(112) 

which is obtained by inverting its Stieltjes transform. 

Other formulae for densities of Dirichlet means may be found in Regazzini, 
Guglielmi and DiNunno [42]. 

6) Formula (105) is obtained from (104) by integrating between and oo: 



oo 



1-/^ h.ii/G){x)dx = J- ^*-ie-*^('°s«)ii;(e-^'(«))dx 



-tE{logG) 



E 



1 



(A(G))' 



r(i), hence : 



^ (A(G))* ^ 
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The interest of this formula (113) is the following: it allows, in a situation 
where the law of G is not known but when one knows the laws of Dt{G) and 
Dt{l/G) to show that i?(|logG|) < oo as soon as £;((A(G))"*) < oo and 

E{{Dt{l/G))~') <oo. 

Formula (104), once multiplied by x'^ , and integrated between and oo, leads 

to: 

The formula is also true even when the expectations which appear in this ex- 
pression are infinite. 

7) The resuh: 

(ri(G))*^^^C/ (115) 

with U uniform on [0, 1] (see point 3 iii) of Proposition 1.7) is a consequence of 
(104). Indeed: 

P{{Tt{G)y <e-^) = P(r,(G)<e-f) 

fTt{G){x)dx 
ptBOogG) re-T: 



m 

(from (104) applied when replacing G by 1/G): 



1 



dx = 



t^or{t)Jo r(t + i) no 

We note that, from (64), the family of the laws of Dt(l/G), i < 1 is tight as 
soon as E{G) < oo, since E {Dt (^)) = E{G). 



2.2. Proof of Theorem 2.1 

2. 2. a Point 1 is trivial. We show point 2: 
Since, from FruUani's formula: 

= logA-ii;(logG)+i?(log(^l + ^) 

we have, by changing A in 1/A: 

-^g(1/A) + ^i/g(A) = -log A + £;(log G) 



L.F. James, B. Roynette and M. Yor/GGC and Dirichlet Means 



376 



2.2.b Wc now show point 3 of Theorem 2.1: 

i) From formula (101), after multiplying by t and exponentiating, we obtain: 



p-tE(\ogG) 

S(e-^r.(i/G)) ^ f S(e-xrUG)) 

A* 



(116) 



Then, taking the Laplace transform of both sides of (116) in the variable A, we 
obtain: 



g-tBOogG)^ 



2e 



' -/3A-irt(G)^ 

A* 



(117) 



where denotes the Bessel-MeDonald function with index u and where we 
have used formula 5.10.25 in [.H4]. p. 108 and 119. 

We now use (117) to compute /r^ (i/g) by inverting its Stieltjes transform [53]. 

1 



2m vloV \-u-iT] + Tt{llG) -u + ii] + Tt{llG) 



f r,(i/G))} 



-tB(logG) 



rt(G) 



-e-'^Ki_t{2^uTt{G)e 



e^^K^^t{2^JuTt{G)e-'f) 

(118) 



However, it is well known that (see [34], p. 108 and 109): 
• TT < Arg z < I =^ A'^(z) : 



• I < Argz < TT ^ A'^(z) 



i/r(ze-^) 



where the Hankel functions hI^'' and iJ^^'' satisfy ([34], p. 108) 

+ = 2 J, 
Plugging these relations into (118), we obtain: 

'^t{G)\^ 



/r.(i/G)(-)=e-*^('°««)i?{(i^) ' J,_i(2V^If;M)} 
iij We now prove (104), which is a consequence of (102). Indeed, from (102) 



and (103): 



/rt(i/G)(w) - e ^ ^U~5r~ 



E 

fc=0 



T{k + l)T{k + t) 



(wrt(G))'+ ^ } 



g-t£;(logG)yt- 



OO 

ME 



(-1)' 



fc=0 



r(fc + i)r(/c + t) 



(«7tA(G))'} 
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from (53). But, since E{'y^) = ^-j^^^, we have: 



\ ' 1 n 



e 



fc=o 

-t£;(logG) 

r(t) ^ ^ ^ ^ 

On the other hand, as we akeady noticed in point 6 of Remark 2.2, formula 
(105) may be obtained by integrating from to oo in (119). 

2.2.C We now prove point 4 of Theorem 2.1. 

i) Formula (106) is immediate. Indeed, since Tt{G) = ^tDt{G) we get: 

.,e-.-,=„.,(-W.(A,)=.(^^^^) 

Replacing then A by 1/A in this latter formula, we obtain: 
1 

E 



(A + A(G))^ 



= exp{-t(V'G(l/A)+logA)} 

= exp{-t(V'i/G(A)~i?(logG))} (120) 

from (101). We note that (105) may also be obtained by taking A = in (120). 

a) We now show (107). For this purpose, we shall prove that the two members 
of (107) admit the same Sticltjes transform with index t. Indeed: 

after making the change of variables y ~ 1/x 

^g-tB(iogG)^(g-Ar,(G)) from (56) 

= exp-i(£;(logG)+V^G(A)) (121) 

whereas: 

1 



fD,il/G)ix)dx = — / —t fD,{l/G)ix)dx 



(X + xY '-^'^'->' ' A* Jo (1 + f)* 

= exp{-i(logA + ^i/G(l/A))} (from (56)) 
= exp{-t(S(logG) + VG(A))} (from (101)) 

(122) 

The comparison of (121) and (122) and the injectivity of the Stieltjes transform 
of index t imply (107). 
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in) Wc now show (107), for i = 1. by using (110). For t — 1, (110) writes: 

,^,(,) ^ M-Fi/oi-)) ^-^((,og|.-i|)i.,,,.) (123) 

TT 

Thus: 

= cxp|-i?(^log(.G|---|j.WG 

= 1 p--E(logG) 



g-^(iogc)_^^^^^^(l/^)^ from (123) 



since sin(7rFG(x)) = sin(7r(l — FG(a;)) = sin(7rFi/(3(l/x)) . It is possible, by 
using (110), to extend this proof for all t < 1. We leave the details to the 
interested reader. 



2.3. A complement to the duality theorem 

Here again, G denotes a r.v. such that i?(|logG|) < oo, and we recall that, for 
any t < 1 Yt denotes a Bernoulli r.v. with parameter t (see (41) and (42), with 
m = t and m' — Vj . 

Theorem 2.3. For all t e [0,1[: 

1) The density /^j^^^-) may be expressed in terms of G, as: 

sm{TrtFG{l/x)) 
4,(^)(-) = 



a;>0, (124) 



and when 1/x is in the support of Fq. Otherwise replace sin(7ri i<G(l/a;)) by 
sin(7r (1 — t)) when x > 0. 

2) The following duality formula holds: 

f (o^{-)-^^^^^^^^^4^■-'-'f r o (-) (125) 

' sin(7rt(l-FG(l/a;))) Di[jiy^) \x ) ^ ' 

3) Let At : [0,oo[ — > [0, 1] be defined by: 

, ^ 1 / sin (TTt) \ ,^ 

At (w := 1 arctg , \ ^ 126 

^■^^ TTi '^Vcos(7ri) +y/ ^ ' 
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Then. 

( ^oA\ 



FoiM^) = A, ( ) (127) 



4) Equivalently, 



, , , ( E\{Dt{G) - x)-.^ 



\E[{x ~ Dt{G))-, 
Remark 2.4. 

1) We note that the right-hand side of (127) depends on t (t e [0, 1[), whereas 
the left-hand side does not depend on t. 

2) Since (see (43)) Di{^) • A(l/G), the knowledge of the law of 

Dt{l/G) and of that of Dt{G), for one t < 1, allow to determine that of G. 
We shall exploit this fact, in points 2.5 and 2.6 below, to determine the Thorin 
measure of a Pareto distribution and of a power of a gamma variable. 

3) We note that finding an explicit Thorin measure of an arbitrary GGC is 
akin to finding the Levy measure of some infinitely divisible random variable. 
Bondesson([ ], Theorem 4.3.2, p. 61), using inversion techniques, obtains an 
expression for the Thorin measure, but notes that it seldom yields explicit ex- 
pressions. On the other hand the use of statement 3) of Theorem 2.3 will often 
lead to tractable expressions for the Thorin measure. 

4) We shall prove (see Section 3.1.b) that At is the cumulative distribution func- 
tion of a r.v. Zt which we shall describe. On the other hand, some trigonometric 
computations allow to see that A^^, the inverse of At in the sense of composition 
of functions, equals: 

1 / ^ sm(TTtx) r 1 , ^ 

= ■ / ,n ^\ ^ ^ 1] (128) 
sm(7rt(l — X)) 

5) Point 1 of Theorem 2.3 is due to James (see[28]). 

2.4. Proof of Theorem 2.3 
2. 4. a We first prove point 1. 

As this point has already been established by James (see [28]), we shall only 
give a quick proof. 
We have, from (56): 



E 



1 



exp(-i^G(A)) = ii?(^^^) 



Hence, changing A in j and using (101) 
1 



E 



exp{-tVG(l/A)-logA} 
exp((t-l)logA) •e*^(i°s)exp(-<Vi/G(A)) (129) 
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We then compute / / ^ \ by inverting its Stieltjes transform: 



/ 



e*s(iogG) 



im|cxp[((t- l)log(-u- i?7)) - tipi/a{-u - ir]) 



2i7r j)io 

-exp[{it - l)log{-u + tri)) -ti^i/ai-u + tv)]} (" > 0) (130) 
It then sufRces to observe that: 

exp {t - 1) log {-u - — > cxpj {t - 1) log \u\ - iTr{t - 1)} = u^-^e.-"^*-'^^ 

and: 



lim exp{(t - 1) log {-u + iri)} = u*-^eJ'^^*-^^ 



as well as: 



t'lpi/ci-u-iv) =tE{\og{l-uG-iriG)) 
= tE{\og{l-uG~iijG)luG<i+log{l^uG-irjG)luG>i) 
^ tE{log{\l~uG\)luG^i^nrP{uG>l)) 



whereas: 

t^Pi/G^u + i7j) ^ t E{\og{\l - uG\)luG^i + iT^PiuG > 1) 

Then, plugging the values of these different limits in (130), we obtain point 1 of 
Theorem 2.3. 

2.4.b We now prove point 2 of Theorem 2.3. 



We deduce from (104) that: 
,)(-) 
sin(7rti^G(l/a;)) 



sm{TrtFG{l/x)) 



x*-iexp{ - t ij(log(-|l - xG\) l,G#i) } 
x*-ie*^('°s«'cxp{ -fii;((log|l-xG|)l,G^i)} (131) 



whereas: 



Di 



.exp{-t^;(log(|i-l|)w,)} 
^ smintFyaix)) ^t-i^tEiiogG)^^^^ - i£;(log|l - xG\)1,g^,)} (132) 
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hence, point 2 of Theorem 2.3 follows, by comparison of (132) and (131). 

2.4.C We now prove point 3 of Theorem 2.3. 
From (125), we obtain: 

4i(^)(^^ _ 8m{TrtFG{l/x)) 



Ai'iFcil/x)) 



with Ar' (x) = '^''^7^''^ , (see (128)) 
sm7ri(l — X) ' 



Then, inverting this formula, since At o A^ ^(x) = x, x [0, 1] we obtain: 



2.4.d We now prove point 4 of Theorem 2.3. This result follows by using the 
generic form of the density of a Pt,i-t random variable multiplied by an indepen- 
dent random variable in the particular cases of (3t,i-tF>t{G) and Pt,i-tF>t{l/G). 
The result is completed by applying the identity in (107). 



2.5. Computation of the Thorin measure of a Pareto r.v. 

Here is a first application of Theorem 2.3. 
2. 5. a Let to > fixed and: 

Xe:=— (0 < 6* < 1) (133) 

7m 



with density: 



- ^(^ + "^) x'-'E{e-^'">+-) (135) 



r{e) ■ r{m) 

The r.v. Xe is a GGC r.v. (see Bondesson [^>], p. 59 ). The rationale of our work 
is now the following: 

• We first compute the Thorin measure associated with Xe- 

• Then, letting 8 converge 1, we shall obtain - as the Thorin measure depends 
continuously (for the narrow topology) on the law of Xg (see Bondesson, [5]) - 
the Thorin measure associated with the r.v. Xi = i.e. to a Pareto r.v. with 
parameter m{m > 0). 
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2.5.b Thorin measure associated with Xg, 9 < 1. 
Since, from (134): 



sup-^ a > ; lim 



a;— a;"~i J 

we deduce from (30) the existence of a r.v. Gg, such that £'(log''' (^)) < oo 
and such that Xe is a {9, Ge) GGC r.v.: 

Xe ^'^^'^ TeiGe) ^'^^^ jeDe{Ge) ^'^^'^ ^ (136) 

7m 

hence: 

DoiGef^^'— (137) 

7m 

On the other hand, from (104): 



Comparing (137) and (135) yields: 



JE(\ogGg 



T{9 + m) 

r(m) 



(139) 



Do{l/Gg) ^'=' 7e+m (140) 

Since, on the other hand: 

dJ-^) ''-'pe,.-e-De{l/Ge) '"^^ Po,i-o ■ 7o+m (141) 

' Gq \ {law) . , [law] 1 

7r7 

we easily deduce from these formulae that: 



V Ffl / 7m, 



The Thorin measure of Xq, which equals: 9PQg{dx), where Pcg is the law of 
Go, is then, by applying (127): 



oo ,,m+S — 1 

dy 



-yVl 



Fg, [-) = Ao ( jl^ I (145) 







where, to obtain (145), we have used (139). 
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2.5.C Thorin measure of (m > 0). 

By continuity, the r.v. ^ is a (1,G) GGC r.v. and its Thorin measure is the 
law of G whose cumulative distribution function is obtained by letting 9 tend 
to 1 in (145). To obtain this limit, we shall develop several computations. 



2.5.i) Development of J. 



(V-Z)" 



dy as /3 = 1 — I 







y 



m+e-i 



e 3 du [y = z + u) 



lo 

-z poc 



= -— / e-"M'^[(z + u)™"'^- (m-/3)(z + w)™"^~i]dM 
P Jo 

with (/3 — 1 — 0) and after integrating by parts: 



u 



z + u 



l-(m-/?)^) 
z + u/ 



du 



V z + u/ 



du 



/3 Jo 

+6-/ e-^iz + ur(\og-^ + ^—)du + o{(3) 

1 



' (l + -) "Ylog-^ + ^) du + 

V z/ V 1 + u z + u/ 



P JO 

since e~'^{z + — j^^du — by integration by parts. Finally: 



dy 



[y - zf " (3^0 P 



+ e-^z™a„(z) + o(/?) (146) 



with 



Cm{z) = I e 
lo 



u \ "V, u 1 1 , 

1 + - log-— + )du 

z / \ I + u z + uy 



2.5. a) Development of z" e-^^x^+'^^^l - lO^^'dx as /3 1 - 61 — >0 



(147) 



/ e-''^x"'-'^{l-xf-^dx 



3™(l-u)'""'5u''du (a; = l-u) 

e^"((m -/?)(!- u)'"-''-^ - 2(1 - u)'"-'3)u'^du 
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with (3 — I ~ 9^ and after integrating by parts 

u \P f m — P 



du 



1^ Jo \l-u J Jo 



z log 



1 — u \1 — u I \ — u 

-z ^ra 



du + o{P) 



with 



— — + e--z'"C„(z) + o(/3) (148) 
P 

[\-"{l-ur(--^+(-^-z)\og-^)du (149) 
Jo V 1 — u VI — u / 1 — u/ 

since: e^"(l — u)'" (y3j; — z^du = 1, by integrating by parts. 
2.5. Hi) Let 

Qi-pi^) ■■= yr (150) 

/ e-'''x'"+''-\l-xfdx 
Jo 

Thus, we have: 

l^z"^ + e-''z^Cm{z)+o{P) 

= (151) 

-j-+e-^z^CUz)+o{(i) 

where Cm{z) and Cm{z) are given by (147) and (149). Hence: 

qi.fiiz) = 1 + /3(C„,(z) - C^{z)) + o{P) (152) 
Plugging this expression in (145), we obtain: 

1 / sin(7r 

= 1 arc tg ' 



{7t9) + 1 + P{C^{z) - Crniz) + 0(/?)) 

1 / sin7r/3 

arc tg 



^(1 - /3) V-cos7r/3 + 1 + /3(a„(z) - C,n{z) + o(/3)) / 
l--arctg( ^ Wf(1/z) (153) 



where F is the cumulative distribution function of the Thorin measure associated 
to the Pareto r.v. with parameter m > 0. 
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1 /a 

2.6. Thorin measure associated with , < a < 1 

2. 6. a Let Xa ■ '''='''' 7^" with < a < 1, with density: 

fxjx) = ax^^ie""" = ax''-^E{e-^^-) (154) 
where 5*^ is a positive stable r.v. with index a: 

i;(e^^'^°) =exp(-A") (A > 0) (155) 
From Bondcsson ([l], p. 60), Xa is a GGC r.v. Since: 

a = supli/ > 0, Km Z^fiM = qX 
Xa is a (a, Gq) GGC r.v., for a r.v. Ga such that: £'(log^(l/GQ)) < oo. Thus: 

Xa ^="'^71)'/" = ro(G„) = 7a^a(Ga) (156) 

We shah now devote some effort to finding the law of Ga, making use in partic- 
ular of formula (107) 
From (104): 

^ ' 1 (a) 

we deduce, by comparison with (154): 

e"'^'-^°<^^''^ = aT{a) ^T{a + 1) (158) 

and 

Da{l/Ga)^'=^ Sa (159) 

Since: ^ 

-Pi ( „ ^ ) = /3Q,l-a • -Dq,(1/Gq,) = /3q,i_q • S'a 

we deduce that the density of ( g V ) equals: 

= y"-^ U-)dx (160) 

where /a denotes the density of Sa- On the other hand, from Chaumont-Yor 
([9], p. 112): 

(161) 

Oa 

hence, since: 

7l (law) (law) (law) (''^a\ /i Gr,\ 

= la-Da(Ga) = Pa^l-a ' llDa[Ga) = 7l ' -Dl ( ) (162) 
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we get: 

(163) 
(164) 

Then, applying (127), we get: 

where, to obtain (165), we have used (157). This expression (165), provides us 
then with the expUcit form of the Thorin measure of (71)"'^^", which is equal to 

2.6.b. We shall now give a more suitable expression of (165). From the relations 
(162) and (163), we deduce: 

P^,,_^.D^{Gy^^'^DJ^)^l- (166) 

hence: 

-,.fJ-) = '^y-' r fo.ia.^i^)dz (167) 

r ^y^ Jy [z-y) 

But, the density of Da{Ga) niay be computed from that of ( 1 /Gq ) , thanks 
to (107): 

fD^iG^){x) = e"^('°s«°)/p„(G„)(l/x). 
Now from (159) and (158), we obtain: 

fD^(G^){x) = r(a + l)/„(l/x). (168) 

We now note the interesting relationship concerning the law of Sa- 
Lemma 2.4 

Let < a < 1, and S denote a positive random variable with density {f{y),y > 
0) such that: 

yfiy) = cE[^^^Ms<y}] (169) 

for some C > 0. Then, S is a stable (a) variable; precisely: 

i;(exp{-Ag}) =cxp{~ ^^^^^~"^ A"}. 

We postpone the proof of the Lemma for the moment, and we note that, 
plugging (169) into (165), we obtain: 
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= A„ 



1 

= 1 arc tg 

7ra 



sin(7ra) • E{ , Jg,^ ls„<y) 



, (170) 



_cos(™) • -B( (^^i-^ ls„<„) + E{^^-^ ls^:,y) 
from (126). 

Proof of Lemma 2.4 From (169), we take the Laplace transform of both sides: 



/ 

Jo 



(y - sr 



dy 



= CE 



dy 



= Cr(l-a)A"-i£;(e"^^). 



Denotmg 0(A) = E{e ^^), we get 



from which we deduce: 



b'{\) =Cr(l-a)A"-V(A) 



0(A)=exp{-^I^ii^A"}. 



3. Explicit examples of GGC variables associated with the 
(Gc, < a < 1) family 

All the examples discussed in this Section are related to the r.v.'s. (Gq, < 
a < 1) introduced in [4]. Below, we indicate the properties of these r.v.'s which 
we shall use. We also recall our notation: 

rt(G„)=7t-A(G«) and (171) 
£;(e-^r*(«°))=exp{-i^°°(l-e-^^)^^(e-"«°)} (A,i>0) (172) 

3.1. The family (Ga, < a < 1) 
(see [4]) 



3.1. a. For < a < 1, the density /g„ of Gq equals: 
asin(7ra) ^"-^(l - a;)"-^ 



/g„ (2;) 



(1 - a) TT (1 - xf" - 2(1 - xYx°' cos(7ra) + x 



2a 



l[o,i](^) (173) 
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In particular: 

• for a — 1/2, Gi/2 follows the arc sine law: 

/gi/2(2;) = - / l[o.i](a;), Gi/2^'^Vi/2, 1/2 (174) 

• for a = 1, Gi is uniform on [0, 1] (175) 

. for a = 0, Go ' t— (176) 
1 + exp(7r6) 

where C is a standard Cauchy variable. 
In general, for < a < 1 one has: 

= (1 + A)" - A" (177) 
The density /ri_„(Gc) of ri_a(GQ) equals: 



a 1 

fr,.^{G^)(.x) = p(Y^ ^ " ^ l[o,oo](a:) 
which may be translated as the following identities in law: 



(178) 



, (7a!i.) 7i_a (iau.) 7i_Q 

ri-o(G«) - 73^7 - ^ (179) 

S.l.b. We note, for < /i < 1, S"^ and 5*^ two independent copies of positive 
stable (/z) r.v.'s, i.e.: 

E{e-^^'^) = exp(-A^) (A > 0) (180) 

and, we let: 

Then (see ['■]:]] or [S], p. 116), the density /^^^ of equals: 

, , sinfTT/i) 1 , , , , 

= -ijr .^ + 2.cos(.,) + i (1^2) 



and we have: 

(!ou.) (Zi_q)-^/° 



1 + (Zi_„)i/a 

or equivalently: 



(183) 



1 {law) 1 {law) /ToA\ 
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We note that the cumulative distribution function Fz of Zn equals 



1 



TT/i 



sin(7r/i) 



Fz^ {x) = 1 arctg ^ /^^ x>0 (185) 



cos(7r/i) + X 



and that its inverse, in the sense of composition of functions, is given by: 

F^\x)^ Bin(.^x) (o<,<i) (186) 

^ ^ sin(7r/i(l -x)) ^ - - ^ ^ ' 

(see Remark 2.4, point 3). 

3. I.e. Although this will not be used in the sequel, we indicate a realization of 
the r.v. ri_Q,(GQ,) which has been at the start of [4]. Let (i?t, t > 0) denote a 
Bessel process starting from 0, with dimension rf = 2(1 — a), with < < 2, or 
equivalently < a < 1. Let, for any t > 0: 



5^"^ := sup{s <t; Rs^ 0}, d^"^ := inf{s >t, Rs = 0} (187) 

'law) 

= 7i an cxpi 
from {Rt, t> 0). Then: 



and let = ^71 an exponentially distributed r.v., with mean 1, independent 



ri_.(G.)''=^4"^-ff^) (188) 

A more general study of quantities such as the RHS of (188), has been developed 
by M. Winkel [54]. 

3.2. Study of the subordinators (rt(Gi/2), t > 0^ and 

(rt(i/Gi/2), t > 0) 

In this Section, G1/2 is a beta (1/2,1/2) variable, i.e.: its law is the arc sine 
distribution, (see (174)). 

3.2.a. Theorem 3.1. 

Let (rt(Gi/2), ^ > 0) denote the subordinator characterized by: 

^(e-Ar,(Gi/.)) ^ g^p| _ i ^°°(i _ g-A.) ^ S(e-"«V2)| (189) 

The following explicit formulae hold: 
1. Laplace transform of Tt{Gri/2)- 

i^(e-Ar,(Gv.))^exp{-t^"(l-e-^-) ^/o(|)e-f} (190) 
= {VT+\-Vxf = = (l + 2A + 2v/A(l + A))^* (191) 
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where, in (190), Iq denotes the modified Bessel function with index 0. [see [34], 
p. 108) 

2. Laws of rf(G-i/2) and of -Dt(G-i/2)- 



(law) 7t n \ 1 C'^"') 1 , 7t+l/2 

= 75 ' M(<bri/2) = 75 = IH 

Pl/2,l/2+t Pl/2,l/2+t 7l/2 



The density of equals: 

3. Wiener-Gamma representation of rt(Gi/2)- 
for any t > 0: 

Here is the dual version of Theorem 3.1: 

Theorem 3.1*. (CifareUi and MeliUi, [')]) Let (rt(l/Gi/2), t > O) denote the 
subordinator characterized by: 

^(e-Ar,(i/Gi/.)) ^ ^^p| ~t - e-^-)^ £;(e"^)} (195) 

Then: 

1. Laplace transform of rt(l/Gi/2)- 

^(e-Ar.(i/G,/.)) ^ (_J__)'* (^^i > 0) (196) 

2. Laws of rt(l/(Gi/2) and of A(l/Gi/2). 

t+l/2,t+l/2: A(l/Gl/2) = A+1/2 ,t+i/2 (197) 



r(i/2 + t)0F 

3. Wiener-Gamma representation of r((l/Gi/2)- 
For any t >Q: 
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Remark 3.2. 

1) Theorem 3.1*. has been obtamed by CifarelH and MeliUi [0]. It would be 
possible to prove Theorem 3.1 by first using Theorem 3.1* and then by applying 
the duality Theorem 2.1. In fact, we shall operate conversely, as we shall first 
prove Theorem 3.1, then we shall show that Theorem 3.1* may be deduced from 
it, due to the duality Theorem 2.1. 

2) By comparing formula (190) with formula (76), we deduce: 

iJr,t>of-' (ir,(G,/,),t>o) 

3) We come back to (196) and we write: 
9 1 /■oo 



i + VTTx i + i(VTTA-i) 

= i?(exp-A5())(i)) (200) 

where (^s[^^\(t), t > O) denotes the stable (1/2) subordinator, Esscher trans- 
formed, with the Esscher transformation with parameter 1 (see [4-5]). In (200), 
e denotes a standard exponential r.v. independent from (^s[^^]^{t), t > O). Hence: 

)"=i.(cxp(-A5«(f))) (201) 



1 + v^TTA 



where the subordinators (^s[^^^^{t), t > O) and (72*, t > 0) featured in (201) are 
being assumed independent. This formula led James and Yor to consider, more 
generally, the subordinator (^s[''^\{'jt) , t > O). In the article [31], the following 
is obtained: 

(rt(Gi/2), t > 0) ''^'^ (5i/2(72t) + S[)li^2t), t > 0) (202) 

where, on the right-hand side of (202), the three subordinators S1/2, s[^^\ and 7 
are assumed independent, and (5'i/2(t), t > 0) ^'^^ (S'i/2(0j t > Oy 
3.2.b. Proof of Theorem 3.1. 



i) We already prove (191). In fact, (191) follows immediately from (176), 
since: 

£.jg-Art(Gi/2)] = (£;[g-Ari/2(Gi/2)])2t 

= (VTTa-Va)''= f— ^ — ^V* 



l + 2A + 2^A(l-f A) 
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a) Wc now prove (190), which is equivalent to: 

£;(e--«V2) = e-*/o(|), ^> (203) 
From the Lipschitz-Hankel formula (see [52] or [38], Th. l.l): 

ly / e-""/,(a;) — = (a + Va^ - l) " (a > 1, > 0) (204) 
Jo ^ 

we deduce: 

\ dx 



yj (1 - e-^")e^t/^ (|) ^ ^ i _ (2A + 1 + 2yJ\{l + A)) " (205) 
then, letting ^ in (205): 

I (1 - e-^")e-*/o( J ) — = log(2A + 1 + 2VA(1 + A)) (206) 



Hence: 



£;(e-Ari(Gi/.)) ^ exp{- /""(l-e-^")^i?(e-^«V2)| 



1 



VTTx + ^/x 
1 



2A + l + 2v/A(l + A) 



This formula, which is of interest by itself, shall not be used in the sequel of this 
proof. 

Hi) We now prove point 2 of Theorem 3.1. 



For this purpose, we shall use the following property of hypergeometric func- 
tions: (see [34], p. 238). Let, for a, (5, 7 reals (cf [34], p. 239): 

F{a,S,j; z) = —^P^—— f\^-\l - t)^-^-\l ~ tz^dt (207) 



(7><5>0, [z[<l) 

- E{{1 - zps,^-5r") (208) 
Lemma 3.3. Let a, 6, c three positive reals. Then: 



(209) 
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(A>0) 
(210) 



Since E{e 



) 



we have: 



E 



fh.. 



[l-xf-^dx 



B{b,c) Jo (x + X)- 
(by definition of the beta (5, c) law) 

B{a + b,c) 1 1 /"^ - 



(1+fr 



dx 



E l 



B{b,c) B{a + b,c) A" 
T{a + b)T{b + c) 1 
T{b)T{a + b + c) A^ 

hence point 1 of Lemma 3.3, from (208). Point 2 of Lemma 3.3 may be obtained 
from similar argmnents. 

We now end the proof of point 2 of Theorem 3.1. From Lemma 3.3, we have: 



E{e 



r(t + i)r(t + i) 11 1 
r(i)r(2t + i) A^^^^'^+2'^^ + ^^-A 

r(^+i)r(t+i) 1 ^i+yr+i 

r(i)r(2t + i) A* 



(211) 



i + \/r 



since (see [34], p. 259): 

F(a, a + - , 2a + 1 ; z 
Using Legendre's duplication formula ([>)4], p. 4): 

r(^t + i)r(f + 1) = 2"2t^(2t + 1) (^withF 



we obtain: 

£;(e 

from (191). Hence: 



f>i 1 



') = ( 



1 



\/rTA + \/A 

7t ('o«>) 



rt((Gi/2) 



(212) 
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The form of the density of rf(Gi/2) follows easily from (212). 

iv) From Proposition 1.3 and (18), showing point 3 amounts to compute the 
inverse of the cumulative distribution function of Gi/2- Now, we have: 

du 



^ Jo \/u(l 



2 

= — Arc sin {\/x) x >Q 
sothat F^/^Jy) = ^^^'{y) i^<y<'^)- 

3.2.C Proof of Theorem 3.1*. 

i) We prove point 1, in two different ways : 
• A direct proof. 

We deduce, as a particular case of the beta-gamma algebra, that: 

iV^ (law) (law) 

— - 7i/2 - e/?i 1 (213) 

where iV is a centered Gaussian variable, with variance 1, and e a standard 
exponential. Hence: 

2 



Thus: 







(l_e-A-)^£;(e %/2) ^ e( ^ (i ^ e-^^)dx 

X \ Jn X 



— (l-S(exp(-A7.7i/2))) 
" 

(from the change of variable: a; = U71/2) 

^ du /■ 1 
1 







u 



(l + Au)i/2 



= 2(log(Vr+A + l) -log2) 

hence: 



E{e '^-1/2^) = exp|-2ilogf 



2 \2t 



-1 + VT+A/ \2 + \ + 2^/T+\' 
We may also prove (196) with the help of the duality Theorem 2.1. Indeed, 
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from (101): 



cxp -tE{logGi/2) 



A* 



from (191) and since £;(logGi/2) = -log 4, from (104) and (192). Thus: 

£;(e-^^'(^)) 



(2 + A + 2yrTA) Vi + VT+A 

ii) We now prove point 2 of Theorem 3.1*. 
Of course, it would be possible to use point 2 of Lemma 3.3 to make this proof, 
in the same manner that we have used point 1 of this Lemma 3.3 to show point 
2 of Theorem 3.1. In fact, we prefer to use formula (104) of the duality Theorem 
2.1. We have: 

,t~i 



/rt(i/Gi/2)(") 



-tE(logGi/2)^ 



E{e 



iDt{G 



1/2) 



^ u^-^ E(e (from (192)) (214) 



m 

Considering now the Laplace transform of the two sides of (214) we obtain: 

4* 

m 

4* 



2 ' 2 ' 



V(t) 

4*S(t+ i) 



E 



1 



1 + A/3t_^i 



(215) 



= 1 



S(i,i+t) 

Note that, by taking A = in (215), we get: 

4*B(t + ^,t+^) 
S(i,t + i) 

(which may also be recovered from the duplication formula for the Gamma 
function). Hence: 

p(e"^^'(cr7^)) = e\( i )*] ^ £:(e-^'^*ft+i/2.t+i/2) 

LVl + A/3t+i/2,i+l/2^ -I 

The other formulae of Theorem 3.1* are now easily obtained. In particular. 



E{ 



since: F^J-^^^^iu) 



. J,ru , we have, from Proposition 1.3: 
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3.3. Study of the r.v.'s ri_a(Ga) and Ti_a{l/Ga), < a <1 



3.3.a. Theorem 3.4. Let, for < a < I, (TtlGa), t > O) the suhordinator 
characterized by: 



We then have the following explicit formulae: 

1. Laplace transform of rt(Gct). 

S(e-^^'(^°)) = ((1 + A)"-A")^ (A,t >0) 

2. Distributions of Ti^ai^a) and of Z3i_q.((Gq). 

(law) 7l-a (law) 7i-q 



(216) 



ri-Q(Ga) 



(217) 

(218) 
(219) 



{where — is, by definition, aPareto r.v. with parameter aj . The density fY-^_^ 
ofTi-a{Ga) equals: 

/ri_„(G„)(a:) 



r(l - a) x^+- 

3. Wiener-Gamma representation of Tt{Gia)- 
For every t> Q and < a < 1 ; 



- (1-e "^) l[o,oo[(2;) 



(220) 



rt(Ga) 



{law) 



(law) 







f 




Jo 








f 




'0 





t—u \\ \ 1/a 



sin(7r(l-a)(^) ) 

sin(7r(l - a)f) 

sin(7r(l-Q)f) N^/' 
.sm(7r(l-a)(^^)) 



d'fii. 



(221) 



In particular, for t = \ — a: 



(law) 



sin(7rit) 



sin7r(l — a — u) 



l/a-i 



d7u (222) 



We note that, for a = 1/2, formula (222) coincides with (194). The dual version 
of Theorem 3.4 is: 

Theorem 3.4*. Let < a < 1 and let (Tt{l/Ga), t > O) the suhordinator 
characterized by: 

^(g-Ar,(i/G„)) ^ ^^p| _ ^ _ e-^n^ E{e-^)} (223) 
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Then: 

1. Laplace transform oi rt{l/Ga)- 

^(g-Ar.(i/G„)) ^ (■ 1 {l + Xr-^ \ (224) 

2. Laws of ri_a(l/Ga) and of Di_„(l/G„). 

ri_„(l/G„)*="'Si_„-t/, i?i_„(l/G„)'="'V (225) 

(we note that the law of Di_q,(1/Gq,) does not depend on a, and it may be 
compared with (197): = U).The density fri^aii/Gc) ofTi-a{l/Ga) equals: 



(226) 



/r.„„(i/G„)(-) - e-^-^)l[o,oo[(-) 

3. Wiener- Gamma representation of Ttil/Ga)- 
For every t > 0: 

r.(l^) /' (227) 



G„/ Jo ^ ^ / sin(^(l-a)(i^)) \l/" 

V sin(7r(l-Q)f ) / 



sini 7r{l— a) 



3i„(,r(l-a)(i^)) 



1/q 



(228) 



Remark 3.5. 

1) Formula (221) was originally obtained by T. Fujita and M. Yor (see [23]). 

2) We deduce from Proposition 1.6: 

£^(^i-c.(^)) =^(G„) = £([/) = i (from (225)) (229) 
Of course, we may verify directly that: 

m^) = \ 

starting from the formula: 

i;(exp(- AeG„)) = \~ ' (230) 

^ 1-a (1 + A)"-1 ^ ' 

(see [4], formula 1.19), then taking the derivative in A = 0. 

3) In [4], Section 4.2.1, the positive variables Xa.i whose Laplace transforms 
equal: 

i?(e"^^-^) = i ^ (0<«<1) 
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have been introduced. We note that, from (224): 

X^/='r,^^{j-) (231) 

4) Since ri_„(l/G„) Si-a • U and ri„„(G„) ^ the relations (221) 
and (227) may also be written as: 



1 _^ (law) 



JJV" _/„ Vsin(7r(l - a - u)) 



VsmlTTll — a — u))/ 



(law) 

and 



l-a 



V sm(7rit) / 



^ ] f sin(7T(l — a—u 
\ sin(-7r?i) 

(law) / d«(Ai ) 



I A. ( sin(TTu) \ 
\ sin7r(l — a— ii) J 

In particular, for a = i, one finds: 



(232) 



3.3.b. Proof of Theorem 3.4. 

i) Point 1) is an immediate consequence of (176) and point 2) of (177), as 
may be shown after some elementary computations. 

ii) We now prove point 3. From Proposition 1.3, it amounts to compute Fq^ 
and its inverse. However, from (182) and (181), we find, successively: 

2-1/" 



1 + zlLl 



P 



I.e.: 



hence: 



.(^)=Fz,_„((y^)") iO<x<l) (233) 



(FZ^ (a-))'/" 
Fg i^) = u (234) 
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then: 

1 



1 - 



sin7r(l — a) /"^ dy 
7r(l-a) io y2 + 2y(cos7r(l-a)) +1 

sin7r(l — a) 



7r(l — a) 



arc tg 



cos7r(l — a) + X 



(235) 



Ai_Q(a;) with the notation of (126) 



and, from (185): 



sin(7r(l — a)x] 



sin (7r(l — a)(l — x)) 
hence, plugging (236) in (234), we obtain: 



1 



= 1 + 



sm 



(7r(l-a)(l-a;))\i/ 



sin (7r(l — a)x) 



(236) 



(237) 



which, thanks to Proposition 1.3, proves point 3) of Theorem 3.4. 

3.3.C. Proof of Theorem 3.4*. 

We use the duality theorem. From formula (104): 



/ri_<,(G„)(a;) 
which we compare with (220): 



,(l-Q)£;(logGc) 



x-^Q r(l - a) 



(238) 



and we deduce: 



i^(logG„) 



log a 
1-a 



(239) 



Then, using this time (101): 

i;(e-^r,(^)) ^ exp{ - i(^(logG«) + log A + V'g„ (^)) } 



from (239) and (217): 



1 (1 + A)"-1- 



(240) 



which establishes (224). The density of ri_Q,(|T-) may be computed from the 

density formula (104) and the knowledge of the law of £'i_ct(GQ) ^ ~ jji/a i 
given by (219). All the formulae in point 2 of Theorem 3.4* follow easily 
from the explicit expression of /rj_^(_i_), and point 3 of Theorem 3.4* follows 
from: 
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and from Proposition 1.3. 

3.4- Study of the suhordinators (rt(/i, + Go), < > O) and 

(rt(^),t>o) 

We recall the following formula, which was established in [4] (see formulae (181) 
and (179) in [4]): 



exp 



{ - /°°(1 - e-^-)^ i^(e-(^+«»))} = A,M>0 (241) 



This formula was obtained in [4] by "letting a tend to in the study of the 
r.v.'s Gq. Here, Go satisfies (see (175)): 

Go , where C denotes a standard Cauchy r. v. (242) 

1 + exp -kC 

3. 4. a. Theorem 3.6. Let ii > be fixed and let (rt(/^ + Go), t > O) denote the 
subordinator which is characterized by: 

^(g-Ar,(p+Go)) ^ c^p| -tj^(l - e-^n^ i;(e-^(^+'^°))} (243) 

Then: 

1. Laplace transform of Ttip + Go). 

j^^^-xr,i>.+Go)) ^ (244) 

2. Laws of ri(/i + Go) and £'i(/i + Go). 

The densities o/ri(/x + Go) and Di(fi + Go) equal: 

/ri(^+Go)(a^) = , M e""" l[o,oo[(a^) (245) 

log(l + -i) X 



log(l + ^) X 

3. Wiener-Gamma representation of rt(/x + Gq). 
For every t > 0: 



(246) 



r(„ + r) ^'"-^ t i + exp(cotg(^)) 

it(Ai + Go) / — — J— / , (7ru\\\ (247) 

7o 1 + /i(l + exp(cotg(^))) 
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Here is the dual version of Theorem 3.6. 



Theorem 3.6*. Let fj, > and {^t{jj^(^), t > O) denote the subordinator 
characterized by: 



(g-^^'^TiTc^)) = exp{ -t r {I- e-^")^ £;(e"^)} (248) 



Then: 

1. Laplace transform of Tt ( ^4*"^^ ) ' 



i.[e--'(^)] = [ilog(i±A(ii^)]' (249) 



2. Laws ofri(^) and D,{^] 



In particular, for ^ ~ 0: 



+ M^)'"-U + , (250) 



r,,±)"L-'.,c, D,{Mu (251) 



The density ofTi { ^^q^ ) equals: 

/p / 1 = ( f C'^pf ^)dy) Mo, oolix) 

3. Wiener-Gamma representation of Ft ( ) • 

Jbr every /i > and t > 0.' 

1 \ (;a;^) /■* l + /i(l + exp(-cotg ^)) 



(252) 



M + Go/ io 1 + exp( - cotg ^) ' ^ ' 

Remark 3.7. 

1) We may take /.i = in the statement of Theorem 3.6*, since i?(log^(Go)) < 
00, but not in that of Theorem 3.6. 

2) Formula (249), with — 0, was obtained by Bondesson ([ri], Ex. 3.3.1, p. 42) 

and the relation Z?i(|l^) ^ = ^ U has been obtained, independently, by Diaconis 

and Freedman ([12], Ex. 7.4, p. 74) and by Cifarelh and MeliUi ([9], Ex. 2, p. 
1393). 

3) The formula Di{^) ^ = ' U may be obtained by letting a > 0, in (225) and 

formula Ti ( |j- J = z.U is formula (225) with a = 0. 



3.4.b. Proof of Theorem 3.6. 

Point 1 follows immediately from (241), whereas formula (245) is a consequence 
of: 
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log(l + i) log(l + i) + « 



1 



log(l + i) 7a+m-1 







2 /"oo 

e~^dx / exp(— a;u)(iw 



log(l + i) 7o 7a+m-i 



log(l + ^) Jo 



On the other hand, formula (246), which yields the density of foiin+Go) follows 
from: 

loeflH J—") / 1 



log(l + i) \l + XDi{^i + Go) 

1 ft dx I ( ^ . . 1 1 A 

after writing 



log(l + i) J_i_ X 1 + \x \ x{l + Ax) X 1 + A.T 

Then, the last point in Theorem 3.6 follows, with the help of (242), from: 



F^IgSv) - ^+i + cxptg(f(l-2,)) 
^ 1 + M(cxp(cotg(7rt/)) + 1) 
1 + exp cotg(7ry) 

Proof of Theorem 3.6*. 

It may be proven by using the duality theorem. From (104) and (245): 

/r,(,+Go)(0)==e^(i°^('^+««)) ^ 



(254) 



I.e.: 

i?(log(M + Go)) - -log (log (l + ^) ) (255) 
We then apply the duality theorem: 

E{e~^^*^'i^^^) = i;(e-irt(M+Go)) . J_ . g-*s(iog(M+Go)) 



(from (255) and (249)) 



log(l + i 



1 p + A(l + .)-- 
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We now compute the density of Fi ( ^^^^ ) by using (104): 



(from (255) and (246)) 



1 1 



M + y 



e"^ dy (256) 



Formula (250) now follows easily from (256). 
We deduce, from Proposition 1.6, that: 



E 



i?([/ + m) = ^ + M- 



In particular, we have: 

E{Go) = E{U) ^ J 

Finally, the last point of Theorem 3.6* follows from Proposition 1.3 and from 
(254). 

3.5. Study of the subordinators (rt(Gi), t > O) and 

(rt(i/Gi), t > 0) 

In this Section (see (175)), Gi denotes a uniform r.v. on [0, 1]. 

Theorem 3.8. Let (Ft((Gi), i > O) denote the subordinator characterized by: 

dx 



^(g-Ar,(Gi)) ^^^p| ^(g-.Gi)| (257) 



Then: 

1. Laplace transform of rt(Gi). 

i?(e--*(«^))=(^^^^)' X,t>0 (258) 

2. Laws of Fi(Gi) and Di{Gi). 

The expressions of the densities o/Fi(Gi) and Di(Gi) are: 



/r 
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sinf?) , , , ^ 

(260) 



/r>i(Gi)(a;) = — T^TTx l[i,oo[(a;) 

Tr[x — 1)^ ^ 



3. Wiener-Gamma representation of rt(Gi) 
For every t > 0: 



r,iG,f-hf thL^'-\f ^ (261) 



u Jq t — u 

Here is now the dual version of Theorem 3.8, which is due to Diaconis and 
Kemperman [13]. 

Theorem 3.8*. Let {Vt (57) i ^ ^ 0^ denote the subordinator characterized 
by: ^ 

£;(e"^^'<'fr)) =exp| -i^ (1-e-^^)^ £;(e-*)| (262) 

Then: 

1. Laplace transform of rt(l/Gi). 

S(e-^i^'(i/^i)) = (e(l + A)-^)* t,X> (263) 

2. Laws of ri(l/Gi) and Di(l/Gi). 

The densities ofTi(l/Gi) and Di(l/Gri) equal: 

e sin(7rx) 1 

/Z3i(l/Gi)(^) = ^^(13^)1^ ^[O'll^"") (2*^^) 

3. Wiener-Gamma representation of rt(l/Gi). 

For every t > 0: 

m/G^f^^^] f\t-u)d^J'"^'\ fudj^ (266) 
t Jo t Jq 

Remark 3.9. 

1) Theorem 3.8*, and in particular formula (265), are due to P. Diaconis and 
J. Kemperman [13]. Formula (265) may also be found in ([s], Ex. 4.4, p. 98) 
corrected with multiplication by e in that formula. 

2) i) A r.v. Z, which takes values on M, is said to be a Luria-Delbriick r.v. if it 
satisfies: 

E{e-^^) = (A > 0) (267) 

and M. Mohle [39] determined the density fz of this r.v.: 

1 1"^ It 

fz{x) = - / e~'^cos{xt + t\ogt)dt (x € R) 
I" Jo 

= - / e-=^*-*'°st(sin7r<)dt (268) 

71" Jo 
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a) On the other hand, it is proven in ([44], § III, 1.3, p. 1251, with a = k = l) 
that there exists a Wald couple {X, H) , which is infinitely divisible and such 
that: 

• His positive and X and H arc infinitely divisible. 

• E{e-^")- E{e^^) = \ A>0 

.i?(e-^^) = (l + A)-(i+^)e\ 

• E[e^^) = {l + \Y+^e-^ (269) 

Thus: 

£:(e^(-^+i)) = (1 + A)i+^ (270) 
hence, from (258), with t = 1: 

ri((Gi)-(l + X)'="''z (271) 

Hi) We denote hy Z a r.v. defined from Z via the following Esscher transform: 
E{ip{Z)) := E{ip{Z) ■ e-^) {ip Borcl and bounded) 
(we note that E{e^^) — 1, from (267)). Thus, we have: 

i;(e-^^) = S(e-(i+^)^) = (1 + A)i+^ 

so that: 



hence, from (271): 



^^'='-(1 + ^) (272) 



ri((Gi) + Z*'=^Z (273) 



We also note that Letac [35] characterized the law defined by (272). 

3.5.b. Proofs of Theorems 3.8 and 3.8*. 

i) We prove (258), which it suffices to obtain for t ~ 1. We have: 

|>»«(lT4TTx)^i°.(T|l)->o.(in) 



whereas: 



Hence: 



/ (i„e-^^)— £;(e-"^0 = / (1 
Jo 2; Jo 



3-^")(l-e-") — 

X 



= -log^+i") (275) 
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from FruUani's integral. (258) now follows from the comparison of (274) and 
(275). 

ii) We now prove (263): 



'(l_e-^")— £;(e~^) = /°°(l^e-^")— ^^'^'^ 



ux 

I - - I 2 

X J Q X J ^ U 



r ^ r ^ ^ e-""da; (Fubini) 
Ji Jo X 



°° du ^ . .. 
— log 1 + - 



(from FruUani's formula (see [34], p. 6, and (3) above )). 



1 1 



= / log {1+Xw)dw = - ((1+A) log(l+A)-A) 



Thus: 



^(g-AFid/GO) ^ exp{-^"(l-e-^-)^£;(e-^)} 

= exp{-i((l + A)log(l + A)-A)} = - 

^ (14-A)~ 

We may also show (263) with the help of the duality Theorem: 

£,(^-Ar,(i/G.)) ^ E[cxp-\r,iG,)j 



A(l + i)-i 



from (258) and since: £;(logGi) = -£'(71) = -1. 

Hi) Point 3) of Theorem 3.8 and 3.8* follows easily from Proposition 1.3 and 
from the fact that Fg^ (u) = u (0 < u < 1). 

iv) We now prove (264). For this purpose, we use the result of Diaconis and 
Kcmperman [13J, see also [42] , 

, , esinTra; 1 , . 

The relation: 

ri(l/Gi)=7i--Di(l/Gi) 

now implies easily (264). 
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v) We show (259) with the help of formula (104) in the duality Theorem: 

e 

g-xy y \, 277 

TT yy{i-yy~y 



from (276). Finally, in order to show (260), we apply (107): 



1 esin(^) /1\ 

1 — (from (276)) 



X J 1 

— l[oa] 



ex TT \x 

■n[x — 1) ^ 



Appendix 



Interpolation between the subordinators {Tt{^/G), t > O) and 
{.It, t > 0) 

We denote, for every w > 0, by CTu : M — > M the decreasing function defined by: 

, , X sinh u + cosh u x tanh it + 1 , , 

au{x) := — — = (278) 

a; cosh u + smh u x + tanh u 

Since the image of by cr„ is equal to ]tanh u, coth u], then for every positive 
r.v. G, and every m > 0, we have: £'(|log(T„(G')|) < oo. 

Let (rt(au(G)), t > O) denote the subordinator defined by: 

S(exp-Art(a„(G))) - exp{ - (A)} 

= cxp{-tj (l-e-^^) ^S(e-"""(^))} (279) 

Since ao{G) = ^ and cr oo(G) = 1, we have: 

(Ft ((70(G)), t>G) = (rt(l/G), t>0) 
(rt((T^(G)), t>0) = (7t, i>0) 

Thus, the family of subordinators (r((o'u(G)), t > O) interpolates, as u de- 
scribes K+, between (rt(l/G), t > O) and (74, t > 0). 

The aim of this appendix is to show that one may compute "explicitly" the 
Laplace transform, the density, and the Wiener-Gamma representation of the 
r.v.'s rt((T„(G)) in terms of those of (rt(G), t > O), for every u > 0. 
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Theorem A.l. We have, for every u >0: 

1) - ?/>g(o-«(A)) + ^o-„(G)(A) = log (A cosliM + sinhit) + fc^ 

G 



with k,i = E 



log 



G sinh u + cosh u 



2) K'iG)iy) = ^u{Fa\l-y)) (0<2/<l) 



(280) 

(281) 
(282) 



r*K(G))^'= 



(law) 



a^F^^l-l)] Jo a^F^'i 



(283) 



3) The density /rt(o-„(G))of rt((T„(G)), equals: 



/rt(<T„(G))(a;) 

„— x tanh u—t k.u 



El e 



-rt(G)tanh« ( ^t{G) 



cosh M I 

^ — x tanh t 



t-1 



'2v/^r;(Gy^ 



coshi 



r(t) 



E 



(coshu + A(G)sinhu) e .:o=h„(oosh„+Dt(G)sinh„ 

(284) 



Remark A. 2. 

1. For M = 0, we have: cto(G) = ^, = -E'(logG) and the precedmg formulae 
and indeed those of Theorem 2.1. 



2. As u — > (X), the relation (283) becomes, by passage to the limit: 

rt(croo(G)) =" / djs^lt 
Jo 

and it is not difhcult to see, starting from (284). that: 



(285) 



/r*(<r„(G))(a;), 



n— ►oo r(i) 

3. More generally than (280), if we take: 



, , ax + b , dx + b , 

c[x) := and cr(A) := , we have 

cx + a cx + a 

-tPa (5(A)) + i/'^(G)(A) = log (cA + a) + fc(a) 

G 



with fc(cr) := 



aG + 6 



(286) 
(287) 
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We observe that (286) and (287) follow immediately from: 



-V'G(5^(A))+V.(G)(A) = -i? 



log 1 



G 



E 



logfl 



A 



= E 
E 



log 1 



\{cG + d) 
aG + b 



- E 



log 1 



dX + b 1 



cX + a G 



log 



(aG + b) + A(cG + d) (cA + a)G 



aG + b {cX + a)G + dX + b 

{cX + a)G 



E 



log- 



aG + d 



Proof of Theorem A.l. 

Point 1) is a particular case of point 3 of Remark A. 2 and point 2) of Theorem 
A.l is trivial. Let us prove point 3. From (280), after multiplying by and 
exponentiating, we deduce: 



E 



,-Art(a„(G)) 



1 



(A cosh u + sinh u)* 



(288) 



then, multiplying each side of (288) by exp(— /3A) and integrating in A from 
-tanh u to +oo, we obtain: 



E 



tanh u 



„-rt(G)fT„(A)-/3A , 



dX 



tanh u 



(A cosh It + sinh-u)* 



(289) 



Making, on the LHS of (289) the change of variable: A cosh w + sinhu = /i, we 
obtain: 



E 



1 



,(tanh«)(rt(tj„(G))+/3) 



coshu 



E 



i/3 + TtiauiG)) 

, — Pf (G)tanh u-\-l3 tanh u 



dfi 



Thus, after simplifying by cxp(/3 tanh u) and using formula 5.10.25, p. 109 in 
[34], we obtain: 



E 



1 



tanh«rt((j„(G)) 



2e 



-tfcu 



coshii 



E 



(/3 + rt(a„(G)) 

:^t(G)tanhti / ^ 

(G) 



' 2^PTt{G) 
coshw 



We then use the inversion formula of the Stieltjes transform, and thanks to a 
very similar computation to the one made in section 2.2.b, we arrive to the 
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identity: 



„x tanh u f ( \ ^ 771 



coshu 



then to (286), by using the scries development of Jt-i- 
/rt(a„(G))(a;) 

g-2;tanhM-tfc„^t-l xDt(G) >, 

= — E\[cos\\u + Dt{G)smhu) e co,hu(co=h„+Dt(G)ta„h„) k 

I (t) ^ J 

Remark A. 3. Formulae (284) and (285), in the more general situation when: 
(^{G) = fgi^ become: 



1-t 



^ . (^)^ (^)) 

(290) 

where fc((T) is given by (287), and: 

/r.MG))(a:) = ^ e-^^^^'^^'^) i?{(c + dA(G))-* e'^i^mu} (291) 
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List of notations 



G,X,Z denote r.v.'s. 

HG- law of G; fc '■ density of G; Fq : cumulative distribution function of G; 

Fq^ : inverse function of Fq- 

T: GGC r.v. with associated Thorin measure /i. 

m: total mass of fi. 

7a: Gamma r.v. with parameter a (a > 0) ; 71 ^ = ^ e, a standard exponential r.v. 

Pa,b- a beta r.v. with parameters a, 6 (a, 6 > 0). 
U : a uniform r.v. on [0, 1]. 
G: a standard Cauchy r.v. 

Yp-. a BcrnouUi r.v. (0 < p < 1). {P{Yp = 1) = p ; P{Yp ^ 0) = 1 - p) 
Sa- a standard positive stable r.v., with index a (0 < a < 1), and density fa. 
a Pareto r.v. of index S (71 and 75 independent). 

Gs generic notation for random variables considered in sections 2.5 and 2.6. 
(Gq, < a < 1) the family of r.v.'s introduced in 3.1. In particular: 

Gi/2 = Pi 1, Go = — —, Gi = (7 

rm(G) a GGC r.v. with associated Thorin measure m • fiQ. 

{it, t > 0): standard gamma process (i.e.: subordinator with Levy measure: 

{D'r\ <t<m) Dirichlet process with parameter m. 

(Tt{G), t>0): GGC subordinator with Thorin measure ^g, and Levy density: 
•ipG- the Bernstein function of the subordinator (rt(G), < > O): 



^g(A) 



Jo X 



Tt{h) = h{s) d'-fs the Wiener-Gamma integral of h. 
Dt{G) r.v. who; 
Dt{G) satisfies: 



Dt{G) r.v. whose law is characterized by Tt{G) ^'^^74 £'t(G) 



Jo F^ (u/t) 

{Ct, t > 0), {St, t > 0), (Tt, t > 0): hyperbohc subordinators. 

{Jt^\ t > Q), {K!f\ i > 0): subordinators associated to the functions 
and Kjy. 

K^,I^, J^: modified and unmodified Bessel functions with index v. 
S^*\X): Stieltjes transform with index t of the positive r.v. X: 
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(^Sa{t), t > O) standard stable subordinator, with index a (0 < a < 1). 
density of 5„ 5,(1). 
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